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ABSTRACT 
A linear structure can be characterized by its compliance matrix 


C, which is 6x6 symmetrical and positive definite and which relates a 


force 6-vector F, to a displacement 6-vector D, by the relation CF, De 
The inverse, S, of C, is called the stiffness matrix and satisfies SD, = 
F This thesis deals with the problem of finding optimal values of such 
matrices C and S from experimental determinations of a sufficient number 
of vector-pairs (F 5.D,) which are presumed to contain random errors. 

J. E. Brock has introduced this problem area, suggested several dif- 
ferent criteria of optimality, and solved some of the corresponding 
specific problems. This thesis completes the solution to a previously 
unsolved specific problem of this group and contributes computationally 
convenient new solutions to another. Moreover, a computer program, 
originally written for the CDC 1604 has been rewritten, in FORTRAN IV 


Language, as two programs for the IBM System 360 computer, and the 


capability has been significantly augmented. 
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SECTION 1 
INTRODUCTION 
1.1 Source and Statement of the Physical Problem. 
The coordinate system used here is a right handed orthogonal triad, 
whose unit vectors are denoted by ey Ens and e.3 see figure l. 


A generalized force will be denoted by F, where 


eno Pos Man Sas Mang 


is a column matrix. Braces 43 indicate that a column matrix is in- 

tended. The first three elements, Fi» Pos and Fy» are force components 
in the ey eo» and e., direction respectively, and the last three ele- 
ments, Fy» F., and Fes are moment components about these axes. 


Also let a generalized deflection be denoted by D, where 


D = | OP ee De eS we D, 

The first three elements, Di Do» and Da» are displacements in the 
Ey» eos and é., direction respectively, and the last three elements are 
rotation components about these axes, 

It is now presumed that a force F is applied to a (suitably con- 
strained) elastic structure, the application being at point P, and the 
corresponding deflection D is observed at P; see figure l. 

It is postulated that a linear relation holds between the elements 
of F and D, i.e., that there exists a matrix C such that 

CF = D 
It is further postulated that matrix C is symmetric and positive definite, 
and that it has an inverse matrix S (S = owt) which, likewise, is symmetric 


and positive definite. 


Matrix C is called the compliance 
(or flexibility) matrix (of the structure) 
and its elements are called compliance or 
flexibility coefficients. 

Matrix S is called the stiffness matrix, 
and its elements are called stiffness in- 


fluence coefficients. 





Note that matrices C and S depend 
not only on the structure, but also on FIGURE 
the choice of point P. 

There are basically two approaches for obtaining the compliance 
matrix C, or the stiffness matrix S. One employs elastic theory in an 
attempt to come up with the elements of C or S. Even with simple struc- 
tures it may not be easy to write down the elements of these matrices. 
The theory must be stretched beyond certainty and one obtains only an 
approximation to the correct result. 

The second approach, the use of which leads to this thesis problem, 
employs experimental data rather than elastic theory to come out with the 
best C or S. Because of inevitable random errors present in the data, it 
is impossible to obtain the matric C or S, without using some mathemati- 
cal procedure that is able to minimize the influence of these errors, 
and at the same time to assure that matrices C and S are symmetrical and 
positive definite. Furthermore the data must contain a sufficient num- 
ber of independent pairs of vectors F and D. 

The question of how many pairs of vectors are required is consider- 
ed in section 7. However in practice, it is anticipated that many more 


pairs will be available than the theoretical minimum. 


1.2 Brief History of this Problem and Related Problems. 

Professor John E. Brock Whe of the Mechanical Engineering Departe 
ment at the Naval Postgraduate School, introduced this problem. He has 
defined different criteria of optimality and, corresponding to each such 
criterion, has established a matrix equation, the solution of which will 
minimize the corresponding errors and permit determining the desired sym- 
metric and positive definite C or S matrix. 

In the next chapter there will be given the mathematical descriptions 
of four criteria of optimality that Professor Brock has introduced in 
reference [1], and the developments leading to the matrix equations that 
which must be solved for an unknown square matrix. 

A description and an explanation of how to use the computer programs 
which has been developed for solving these equations will be given later. 

The data used to test the programs are not from an actual experi- 
ment. An explanation will be given of how sample data, containing random 
errors, are generated for test purposes. 

In this thesis, matrix C has been restricted to be symmetrical and 
positive definite, as is required in the structural problem under con- 
sideration. Two other related situations will be considered in Appendix 
D. 

In general the different criteria of optimality that will be given 
in the next chapter can be applied to all system$which are described by 
a matrix equation of the form CF = D, where C, F, and D, need only to 


be comformable for matrix operations. 


Eiinie ee in brackets refer to references listed in the Bibliography 
page 43. 


SECTION 2 
MATHEMATICAL DEVELOPMENT OF THE PROBLEM 
2.1 Physical Statement of Criterion One. 


It is assumed that there are n pairs of vector observations F and 


As criterion one, it is desired to determine C so as to minimize 
the sum of the square of all scalar errors, viz., to minimize 


2 


Qe oO (1) 


x, = 


it 
= 


M3 
Ms 


J 
where Lhe = (-th element of vector ff 


Now let 


These are ( Gx ) matrices. 


From equation (1) 


a 


x SE CepBI + (Ce BI. +R] @ 


Expanding the first term of equation (2) 


TL an 2 — 2 = 2 ae » Dp 
v3 Car. it) = ({cF,-D,},) + (icf, -3,},) EEE (ick,-D,, b) (3) 
ahs 
If all the terms in equation (2), are expanded as in equation (3), 
it is easy to see that these terms are nothing other than the diagonal 


terms of the matrix product: (cu-v)! (cu-v), and therefore it is possible 


to write 


10 


AX, = th | (cu-vy (cu-v) J (4) 


where tr stands for the trace of a matrix and superscript T denotes the 


transposed of a matrix. 


2.2 Alternate Criteria. 


D 


F. 
Swe 


FIGURE 2 
As may be seen from Figure 2, in a one dimensional case, criterion 
one requires minimizing the sum of the square of the vertical errors a. 
It is natural, therefore to consider, as criterion two, minimizing the 
sum of the squares of the horizontal errors, b. In the six dimensional 


case that this thesis is concerned with, this means minimizing 


n G eo 2 
KH. => > CLSh-Fs 32) (5) 
jz! isi 


Equation (5) is in the same form as equation (1), with the roles of 


F. and D, interchanged and using the stiffness matrix S instead of its 


a 


Ea 


inverse, the compliance matrix C. 
Using the same approach as before this sum can be written in the 


following way: 


6 rsa [ (sv-u) (Sv- oy (6) 


It is possible too, to minimize the sum of the areas of the "error 
triangles", whose legs, in a one dimensional case, are a. and bis this 


is the same as minimizing 


x, = t [(cu-v)'(u-sv)] 7) 


Now, think also, Of minimizing a weighted sum of the square of the 
horizontal and vertical errors (BX, +/3B2%Xq )- This criterion 


number four is the same as minimizing 


X= th[ A (cu -v) (CU-V) + 3, (sv-u) (s v-u)y} 8) 


where B, and Pr are scalar quantities which may be specified so as to 
designate the relative importance of the two separate parts of which it 
is composed. 


From equation (8), without loss of generality, it is possible to 


make jer = 1 and Pr = jp , and write equation (8) as 
T Tv 
Ay =t[(cu-v) (cu-v) + 6 (SvV-U) (Sv-U) | (9) 


Brock [1] also suggested a fifth criterion of optimality but since 
no progress has been made toward satisfying this criterion, it is not 


included here. 


2.3 Five Theorems on the Trace of a Matrix. 


In deriving the matrix equation corresponding to each criterion of 


V2 


optimality, the following five simple theorems concerning the trace of 


a matrix will be used. The proofs are simple and are omitted. 


2.4 


THEOREM 1: If the matrix P and the matrix Q are comformable in 
either order, then each product is square and 
tr (PQ) = tr (OP) 
THEOREM 2: If the square matrices P and Q may be added, then 
tr (P) + tr (Q) = tr (P + Q) 
THEOREM 3: If P is a square matrix, then 
ac 
tr (P) = tr (P) 
THEOREM 4: If P is a matrix and x is a scalar, then 
T 2 
ae a aw 
aX 2 
THEOREM 5: If P is a square matrix and x is a scalar, then 
Qtr(P) — tr [at 
dx dx 
Note that theorem 1 permits one to write: 


tr (ABCD) = tr (BCDA) = tr (CDAB), and so on. 


Derivation of the Matrix Equations which assure the Optimization. 


In sections 2.1 and 2.2 four different criteria of optimality have 


been given. This section devotes attention to the derivation of a matrix 


equation corresponding to each criterion. Solving any such equation 


assures achieving the corresponding optimization. 


Now let x be any given element of C, and let C* denote 


a) Equation (4) is 


a, =h[ (cu -v) (cu-v)] 


= th (U'CCU-V"CU-UCV+V'v) 


Note that instead of out C has been written since C is symmetric. 


oC 


, then 
3% 


13 


2% =H (U'C*CU+UTOU-VCU-UT'V) 


9x 


= tr [ (cuu™ +UU'C -uv™-vu')( ‘a *) ; 


=t[(cP +PC'-a)(c*)] 


where 


sh 
P = ut, and Q = UV + vut 


are obviously symmetric. 
Thus, if C satisfies the matrix equation 


CP + PC = Q (10) 
ax, 
3x 


which will be minimum, since XK, is the sum of the squares of real 


then = 0, for all elements x of C and there is an extremun, 
numbers. 
Because of the simple relationship between OX, and Xo » 1E oe 
evident that a similar equation applies for criterion two. 
This equation is 
SR + RS = Q (11) 


where 0 is as before and R = wi 


b) Equation (7) is 
x, = te [(cu-v) (U-sv)] 


The solution will be to determine the matrices C which make Xs 
stationary; among these will be one which minimizes X58 
Knowing that S = cnn and both C and S are symmetric, after some 


manipulation 


14 


A, =h (UCTU-V U-UV+V'SV) 


As before, let x be any given element of C, and let C* - 


and S* -23 





Since CS = I, then 
CS* + C*¥S = 0 
and S* = -SC*S. 


Thus 





as = fh (UC *U-V"SC*SV) 


= t[(uu"-svv's)(c*) ] 


This will vanish and Ax will be extremal if S satisfies the matrix 
equation 


uut = svv's 


After post and pre-multiplying by C 
cuu'c = wt (12) 


c) Equation (9) is 


X= tr [(cv-v) (cu-v) if (Sv- u) (sv-u)] 
-h[vV-U'CV-v"CU +UCCU + 


A(vssv-Usv -v'su +uU'v) | 
As before, let x be any given element of C. 


Thus, using notation introduced previously, viz. 


Dole aie en © saul eae 


it may be seen that 


15 


aoe = te Luc*v-vC7U + UTC*CU +UTCOU + 


3B (Vs*sv +v'ss*v - u's*v -V"S*u)] 


fe [ (-vu-uv7+ CULT + uve )(c*) + 
/2 (svv* +vv"s -VUT- UVT)(S*)7 

te [(-Q +CP+PL)(C*)+A(skeRS-Q)¢SC"S)| 

= te [{ CQ+PC+CP)-BS(sR+RS-Q)5 } (C*) J 


iH} 


which will vanish if 
pc+cp-Q=  s(sr+Rs - Q)S (13) 


If the right hand side of equation (13) is set equal to zero, the 
left hand side is precisely the same as equation (10), which corresponds 


to criterion number one. If the left hand side is set equal to zero, 


the right hand side is nothing other than equation (11), which corres- 


ponds to criterion number two. 


16 


SECTION 3 
SOLUTION OF EQUATIONS FOR CRITERIA 1 AND 2 

Equations (10) and (11) that minimize criteria 1 and 2 respectively 
are of the same form, where the matrices C, S, P, R, and Q, are all sym- 
metric. The solution by Brock [1], given here, is for equation (10), but 
is applicable to equation (11), with the roles of P and R interchanged. 
The solution in equation (10) will be for the unknown matrix C, where 
in equation (11) is for the unknown matrix S. 

In the physical problem of primary interest, matrices C, P, and Q 
are all 6 x 6, but, for the sake of generality, it will be presumed that 
they are m x m. 

In equation (10) 

CP + PC = Q 
there are really only N = 4m (m + 1) unknown elements of C. These un- 
knowns can be arranged in a column matrix according to the following 
procedure: 

Use the notation Me =% (j - 1) j. 


Given i, j, let t, . =M + Tintin (educng ye. 


i,] max(i, j) 
For example, with i = 3, and j = 5, then 


t = € 
and tog = t = at) Coy 2 


iS 


A 17 


Table 1 gives values of t. j EO reso ipeieee genera LO 


Next,a column matrix W of N elements is formed from the N unknown 


elements oe of matrix C. Element ay of C becomes the t,. .th element 
> 


of W; note that t = N. 
m,m 


17 





9 | 37 | 38 | 


TABLE 1. Values of ty j 


W = { EL di { CyyreeeSam f 


The column matrix Z of N elements is formed in a similar way from 


the known elements qy j of symmetric matrix Q. 
3 


Z = {2 .2); Zasree styl 


where 9 = 4; /at f Lp 


and Si, is the Kronecker's delta. 

Finally the coefficient matrix P is formed as the sum of m (N x N) 
matrices Poe I = 1,2,...,m, each of order N x N, which themselves are 
formed from the (m x m) matrix P = uu", according to the following pro- 
cedure; 

Po is formed from P by inserting rows and columns of zeros, an oper- 
ation that is denoted by (*)(K)(L) which means insert L rows of zeros 
after the eal row in the original matrix P, and insert L columns of 
zeros after the qh column in the original matrix P. The Ps is obtained 
from P by the operation: 
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(a) (*) (0) (M,) 


(b) (*) (I) (1-1) 
(c) (*) (I+1) (T) 
(d) (*) (I+2) (I+1), etc. 


until such insertions fall outside the range of the N x N size. 


Then: 


and solving the system 


gives the column matrix W, whose elements can then be rearranged to 
produce the desired matrix C. 

Bickley and McNamee [4] have considered a more general equation of 

the form; 

AZ + ZB = F 
and they give four methods of solution, which they have named: The Big 
Matrix, the Irrational Solution, the Semi-Rational Solution, and A Ration- 
al Solution. 

The method given by Larsen [5], and the method given by Brock [1] 
and reproduced here, appear to be additional distinct methods; both, 
however, probably are closely related to what Bickley and McNamee refer 
to as the "big matrix", in which A = B. 

As the results obtained with the method given here is quite satis- 
factory, no further study has been done with the methods given in refer- 
ence [4], but it might be of interest to one working with this problem 


to try suitably specializing Bickley and McNamee's methods. 
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SECTION 4 
SOLUTION OF EQUATION FOR CRITERION 3 
To solve equation (12), postmultiply by vuls then 


cuu'cuut = vw uu! (14) 


cuu' = cw" wut)? 


In equation (14), vu! and wi, both are symmetric, but the product 


(15) 
wi vt is non-symmetric. 

Brock [6], has solved equation (15), using a highly sophisticated 
computer library subroutine, based on a Laguerre-Hessenberg algorithn, 
to find eigenvalues of wwii, which are then used in a method described 


te (wT uv?) ®, 


by Frazer, Duncan, and Collar [3], in order to find CUU 
from which C may easily be obtained. This procedure is available in a 
computer program for the FORTRAN 63 language to be used in the CDC 1604 
computer. 

Owing to the difficulties in converting this program to FORTRAN IV 
language used in the IBM System 360 computer, recently installed at the 
Naval Postgraduate School, the following two other procedures were develop- 
ed and have successfully been applied to solve equation (12). 

a) Iterative procedure for finding the square root of a non- 
symmetric matrix. 

Pulay [7], has given a method that determines the principal square 
root of a positive definite symmetric matrix using an iteration pro- 
cedure. 

Using his idea with some modification, a method for obtaining the 
Square root of a positive definite non-symmetric matrix has been suc- 


cessfully used and applied to several examples. 


In equation (15) 


20 


cou? = (w' wr")? 


4 Ne 


let VV UU A 


and AY =D+B (16) 
where D is a diagonal matrix, whose elements are the positive square 

roots of the diagonal elements of A, and B is a matrix formed from euut o 
D, and @ is an approximation of the matrix C that will minimize the cri- 


terion number 3. 


The matrix © can be obtained from the equation 


CU = V 
éuu? = vot 
and é = wi wut)! 


This is one way Of finding the matrix C, suggested by analogy to 
Pulay's method. However, in the computer program used in this thesis 
to solve equation (15), the matrix € is taken to be the matrix C that 
satisfies criterion number 2. This is done in order to save computer 
time. 

From equation (16) 

ou na 
DB + BD=A - Dp? - Be 

This equation can not be solved for matrix B due to the term B, 
but following Pulay's procedure it can be transformed into a simple 
iteration formula 

ye Sa Dees Be 


and the (i, j)-th element of Bs can be obtained in a simple manner from 


the formula 


(B, a 7D) DD, [(A-D),, ; B.),, | 


aM 


To obtain (Bidig? the iteration procedure must be carried out until 
the difference between two successive approximations of (Bids become 


less than a specified allowable error, i.e. 
B57 By py7 A 
, -10 : .; : 
where A in the actual program is 10 » and the Oi; implies comparing 
element by element. 
After matrix B has been obtained 


C= (D+ B) (uur) 7! 


(wt wy 2 wut) 7} 


OT C 
The procedure described above leads to output 3A in the actual 
computer program. Output 3B is the solution of basically the same equa- 
tion, but instead of equation (12) being postmultiplied by vut, 1G rs 

premultiplied resulting in 


vure = (uv w')? 


and using the same procedure as before the square root of vut wv? can be 
obtained, and therefore the unknown matric C. 

Output 3C is simply the arithmetic average of 3A and 3B. 

b) The second procedure that solves equation (12) uses some 
matrix algebra manipulation, and instead of working with a non-symmetric 
matrix , it requires taking the square roots of two different symmetric 
matrices. This procedure was called to the attention of the writer by 
Professor Brock. 

Observe that starting with equation (12), cuutc = wi and pre- 
multiplying by ut and postmultiplying by U, the following matrix equa- 
tion is obtained 

ut cuu cu = uw y 


and 
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uicu = (u'w"u)? 


where vlw'u is a symmetric matrix, but as U and V have dimensions (6xn), 

where n is the number of pair of vectors available, this formulation of 

the problem requires taking the square root to an (nxn) symmetrix matrix, 

with the disadvantage of dealing with a big matrix, for example (1000x1000). 
Accordingly, return to equation (12) 


cuutc = wi! 


Let uut = A, and without much loss of generality assume that the 
eigenvalues of A are distinct. In this problem A is a (6x6) matrix. 

Then 

A = MLM" 
where M is the modal matrix of A 
and ie is the spectral matrix of A. 

This is true since A is symmetrical, te is simply a diagonal matrix 
of the eigenvalues of A, and M is a square matrix whose columns are the 
corresponding eigenvectors of A. M is orthogonal so that mM: = mt, 

Since ig is a diagonal matrix, it is easy to find its square 
root, L. Theoretically there are a number of different choices of L 
since either the positive or the negative square root of each eigen- 
values can be chosen, but in this application it is convenient and satis- 
factory to use only the positive values. 


sou A =u = Mum! 
mitimt 


EES 


where E = ML is (mxm), i.e., 6 x 6 in the applied problem. 
Let B = wi 


and equation (12) now takes the form 
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cEE‘¢ 


E CEE’ CE = E BE = D 


B 


where D is an (mxm) symmetric matrix. Thus, 
ECE = p? 


Like matrix A, D can be written in terms of its modal and spectral 


matrices 
D = rilerta 
D? = mlm 
Therefore 
E CR = silat 
Now gt = i Ne = Daa 


where ie is simply the diagonal matrix of the reciprocals of the corres- 
ponding elements of L, 
Thus 
c= ue! miata 
This procedure has been written in the FORTRAN IV language to be 


used in the IBM System 360 computer, and in the actual program is output 


3D. 
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SECTION 5 
SOLUTION OF EQUATION FOR CRITERION 4 
The equation (13), that minimize the criterion 4, is nonlinear in 
the unknown matrix C (= gaye and an iteration method is used for its 
solution. 
Given 8 , choose a convenient Pr -f . Now from equation (13), 


consider the iterative scheme indicated by 


CP + PC, --Q =f SMS/R A RSS - Q)S, (17) 
where 
a a 
and C= aa ape? 
et p 
C is the matrix that satisfies criterion 1 
C, is the matrix that satisfies criterion 2 


Now, having os S, is determined, from which the right hand side 


1 
of equation (17) can be evaluated, and the following equation can be 
obtained: 
ae = 0* 
c.P PC, Q (18) 
where 
k = + ata = 
Q Q j 5, (S,R RS, Q)s, 
Equation (18) is of the same form as equation (10), and can be 
solved using the same procedure. 
The C, obtained from equation (18), is used again to obtain S. 


aL 


and a new C 


ar: 


i+1 is obtained using the iterative procedure until two 


successive approximations of C. become less than a specified allowable 


error, viz. 


¢~@). =A 


(C41) 5k - i” jk 


25 


where A in the actual program is taken as Omer and O 5. implies 
comparing element by element. 

When the convergence is achieved, i is increased and the iterative 
procedure is carried out again. This procedure is followed until finally 
convergence is obtained with if = f3 | 

The reason for using values of P initially less than f is to 
assure and optimize convergence of the entire process. For example, in 


one study where (6 = 0.8, it was found convenient to use the sequence 


f =1051j0s20e see One. 
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SECTION 6 
DIGITAL COMPUTER PROGRAMS 
Two main computer programs have been written in FORTRAN IV Language 
to be used with the IBM System 360 computer. These programs are called 
CECI 1 and CECI 2. They make use of several subroutines and to assure 
accuracy the main programs and the subroutines employ double precision 
arithmetic. 


In program CECI 1, the following subroutines are used: 


ADDEM AINVM 
BIGCM BIGSM 
BLANCA JACVAT 
MAUSS3 OUT PU 
PRINM RECIP 
REDEM SQRERT 
SQROOT SQUAM 
TRACM TRANM 
ULTIMA ULTMB 


In the program CECI 2, the following subroutines are used: 


ADDEM AINVM 
BIGCM BIGSM 
BLANCA MAUSS3 
OUTPUT PRINM 
REDEM SQUAM 
TRACM TRANM 
ULTIMA ULTMB 


Both programs make use of a number of functions from the monitor 


library tape which are not listed here. 


ZT 


Program CECI 1 solves the equations that satisfy criteria numbers 
1, 2, and 3. The program gives one output for each of criteria 1 and 
2, and four outputs for criterion number 3. For each of these outputs 
the program prints the compliance matrix C, and for each such matrix, 
computes and prints out three criteria numbers, as follows; 
CRITERION NUMBER ONE  =%n [(cU - V)*(cU - V)] 
CRITERION NUMBER TWO =tr [(sv - ut (sv - v)] 
CRITERION NUMBER THREE = [(cu - v)‘(u - sv)] 
Program CECI 2 solves the equation that satisfies criterion number 
4, prints the compliance matrix C, and gives four criterion numbers; the 
first three are the same as in CECI 1 and the fourth is given by: 


CRITERION NUMBER FOUR =T[(cU - v) (cu ~ V) + B (sv-u)? 


(SV- U)] 

Most of the matrices are regarded as doubly subscripted arrays as 
follows: 

If there are M rows and N columns, the array in the computer memory 
has dimension (M+ 1, N)3 for example: if the matrix has 6 rows and 50 
columns, the dimension of the array is (7 x 50). This extra row permits 
the identification of the matrix for debugging purposes and gives this 
information to those subroutines which use it for processing or monitor- 
ing purposes. 

Let e, refer to the array and let a, . refer to the matrix. Then 


bs | i) 


is an integer identifying the array and the matrix. and e 


“1,2 1,3 


are integers M and N respectively. For the remainder of the elements, 


“151 


es . =a, ,. The other matrices used in the main programs that have 
itl, j i,j 
the dimensions (N,N), are used as scratch matrices. They are needed 


when calling some subroutines that have dummy subscripted variables, due 


28 


to the mode of storage (column-wise) used by the IBM System 360 computer. 
The array for an M-element column matrix has dimension (M + 1), where 
the first element is the integer M. 

The programs can be used in two different ways: (1) for classroom or 
more investigation purposes and (2) for use with data from an actual ex- 
periment. This is done by means of the card called MODE, where MODE = 1, 
2, respectively. 

If MODE = 1, the program will read the matrices C and U from the 
data cards, and generate matrix V = CU. For testing purposes random 
errors are introduced in matrices U and V. The explanation of how these 
errors are introduced will be given in the next section. 

If MODE = 2, the programs will read the matrices U and V from the 
data cards. This data would be obtained from an actual experiment. 

The explanation of how to use these programs, and how the data is 
punched in cards will be given in Appendix A. 

Unless other FORMAT statements are called for, the only printed out- 
put is as follows?- 

a) In main program CECI 1, with MODE = l. 

(1) MATRIX P FOLLOWS 
MATRIX NUMBER 1, 6 ROWS, 6 COLUMNS 
Listing of matrix P in FORMAT (1X, 8D12.4) 
(2) MATRIX BE FOLLOWS 
MATRIX NUMBER 2, 6 ROWS, NN COLUMNS 


Listing of matrix BE in FORMAT (1X, 8D12.4) 


Z 
The nomenclature used in the outputs of the programs CECI 1 and 
CECI 2 is given in Appendix C; see pages 76 and 89. 


29 


(3) 


(4) 
(5) 


(6) 


(7) 
(8) 


(9) 
(10) 


(11) 


MATRIX EB FOLLOWS 

MATRIX NUMBER 3, 6 ROWS, NN COLUMNS 
Listing of Matrix EB in FORMAT (1X, 8D12.4) 
ERROR IN BLANCA = (number given) 

MATRIX U FOLLOWS 

MATRIX NUMBER 4, 6 ROWS, NN COLUMNS 
Listing of matrix U in FORMAT (1X, 8D12.4) 
MATRIX V FOLLOWS 

MATRIX NUMBER 5, 6 ROWS, NN COLUMNS 
Listing of matrix V in FORMAT (1X, 8D12.4) 
ELAPSED TIME = (time given in sec.) 

OUTPUT NO. 1, MINIMIZES V SQUARE ERROR 
DESIRED MATRIX IS 

MATRIX NUMBER 17, 6 ROWS, 6 COLUMNS 
Listing of desired compliance matrix in FORMAT (1X, 8D12.4) 
CRITERION NUMBER 1 EQUALS (value given) 
CRITERION NUMBER 2 EQUALS (value given) 
CRITERION NUMBER 3 EQUALS (value given) 
ELASPED TIME = (time given in sec.) 

OUTPUT NO. 2, MINIMIZES U SQUARE ERROR 
DESIRED MATRIX IS 

MATRIX NUMBER 33, 6 ROWS, 6 COLUMNS 
Listing of desired compliance matrix in FORMAT (1X, 8D12.4) 
CRITERION NUMBER 1 EQUALS (value given) 
CRITERION NUMBER 2 EQUALS (value given) 
CRITERION NUMBER 3 EQUALS (value given) 


ELAPSED TIME = (time given in sec.) 
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(12) OUTPUT NO. 3A, MINIMIZES TRIANGLE ERROR 


(13) 
(14) 


(15) 
(16) 


(17) 


(18) 


(19) 


DESIRED MATRIX IS 

MATRIX NUMBER XXX, 6 ROWS, 6 COLUMNS 

(XXX may have any large integer value, due to the iteration 
procedure used to obtain output 3A) 

Listing of desired compliance matrix in FORMAT (1X, 8D12.4) 

CRITERION NUMBER 1 EQUALS (value given) 

CRITERION NUMBER 2 EQUALS (value given) 

CRITERION NUMBER 3 EQUALS (value given) 

ELASPED TIME = (time given in sec.) 

OUTPUT NO. 3B, MINIMIZES TRIANGLE ERROR 

Same as number (12) 

ELAPSED TIME = (time given in sec.) 

OUTPUT NO. 3C, MINIMIZES TRIANGLE ERROR 

Same as number (12) 

OUTPUT NO. 3D, MINIMIZES TRIANGLE ERROR 

Same as number (12) 

ELAPSED TIME = (time given in sec.) 


TOTAL ELAPSED TIME = (time given in sec.) 


b) In main program CECI 1, with MODE = 2. 


(1) 


(2) 


(3) 


MATRIX U FOLLOWS 

MATRIX NUMBER 1, 6 ROWS, NN COLUMNS 
Listing of matrix U in FORMAT (1X, 8D12.4) 
MATRIX V FOLLOWS 

MATRIX NUMBER 2, 6 ROWS, NN COLUMNS 
Listing of matrix V in FORMAT (1X, 8D12.4) 


ELAPSED TIME = (time given in sec.) 


3.1 


(4) OUTPUT NO. 1, MINIMIZES V SQUARE ERROR 


DESIRED MATRIX IS 


MATRIX NUMBER 14, 


6 ROWS, 6 COLUMNS 


Listing of desired compliance matrix in FORMAT (1X, 8D12.4) 


CRITERION NUMBER 


CRITERION NUMBER 


CRITERION NUMBER 


1 EQUALS (value given) 
2 EQUALS (value given) 


3 EQUALS (value given) 


From here the printed out is the same as when MODE = l. 


c) In main program CECI 2, with MODE = l. 


(1) 


(2) 


(3) 


(4) 
(5) 


(6) 


(7) 
(8) 
(9) 


MATRIX P FOLLOWS 

MATRIX NUMBER 1, 
Listing of matrix 
MATRIX BE FOLLOWS 
MATRIX NUMBER 2, 
Listing of matrix 
MATRIX EB FOLLOWS 
MATRIX NUMBER 3, 
Listing of matrix 
ERROR IN BLANCA = 
MATRIX U FOLLOWS 

MATRIX NUMBER 4, 
Listing of matrix 
MATRIX V FOLLOWS 

MATRIX NUMBER 5, 


Listing of matrix 


6 ROWS, 6 COLUMNS 


P in FORMAT (1X, 8D12.4) 


6 ROWS, NN COLUMNS 


BE in FORMAT (1X, 8D12.4) 


6 ROWS, NN COLUMNS 
EB in FORMAT (1X, 8D12.4) 


(value given) 


6 ROWS, NN COLUMNS 


U in FORMAT (1X, 8D12.4) 


6 ROWS, NN COLUMNS 


V in FORMAT (1X, 8D12.4) 


ELAPSED TIME = (time given in sec.) 


ELAPSED TIME 


(time given in sec.) 


NUMBER OF ITERATIONS = (number given) 
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(10) OUTPUT NO. 4, MINIMIZES V + BETA x U SQUARE ERRORS 


DESIRED MATRIX IS 
MATRIX NUMBER XXX, 6 ROWS, 6 COLUMNS 
Listing of desired compliance matrix in FORMAT (1X, 8D12.4) 
BETA = (number given) 
CRITERION NUMBER 1 EQUALS (value given) 
CRITERION NUMBER 2 EQUALS (value given) 
CRITERION NUMBER 3 EQAULS (value given) 
CRITERION NUMBER 4 EQUALS (value given) 
(11) ELAPSED TIME = (time given in sec.) 
(12) TOTAL ELAPSED TIME = (time given in sec.) 
d) In main program CECI 2, with MODE = 2. 
(1) MATRIX U FOLLOWS 
MATRIX NUMBER 1, 6 ROWS, NN COLUMNS 
Listing of matrix U in FORMAT (1X, 8D12.4) 
(2) MATRIX V FOLLOWS 
MATRIX NUMBER 2, 6 ROWS, NN COLUMNS 
Listing of matrix V in FORMAT (1X, 8D12.4) 
From here the printed out is the same as when MODE = l. 
Actual printouts for both programs, when MODE = 1, are shown in 
Appendix C. 
The listing of the main programs and subroutines will be given in 


Appendices A and B respectively. 
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SECTION 7 
TESTING OF THEORY AND DIGITAL COMPUTER PROGRAMS 

In this investigation the digital computer programs written to solve 
the matrix equations that satisfy the four criteria of optimality have 
not been tested using actual (mechanical) experimental data. The data 
have been generated by assuming a symmetric positive definite matrix C, 
and a set of vector forces {F,} in the form of matrix U. From these two 
matrices the corresponding set of vector deflections { D.} in the form 
of matrix V were obtained as V = CU. 

To simulate experimental data, matrices U* and V* were formed in the 


following way: 


Ux = U + E(U) 


ve = V + E(V) 
where E(U) and E(V) are matrices of errors. 

Then matrices U* and V* were employed as if they were composed of 
experimental data, to be used in the programs to recover an approximation 
matrix C* for C. 

To form error matrices E(U), and E(V), a mathematical procedure was 
devised to perform some alterations in each element of the matrices, U 
and V, in a random manner with an upper limit to the absolute value of 
the percentage of error being introduced. This procedure has not been 
tested for randomness, but the method used to obtain the alterations 


suggests that the values of the elements in matrices E(U) and E(V) are 


indeed pretty good random errors. 


Tee that the use of the asterisk,*, is different in this section 
from its use in section 2. 
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The random errors are generated by a subroutine called BLANCA. 
Because the first row of most of the matrices are used for identifying 
purposes, the subroutine makes use of the element(1,1), which is an 
identifying integer called NUMBER. 

Subroutine BLANCA modifies a given matrix A = [ aij ] by introduc- 


.J, resulting in a modified matrix A* = A+tE. 


ing random errors E = les, 


One specifies the maximum percentage error to be applied. 

First a number Ny = 3.6327 + (matrix identification NUMBER) is 
formed. § is increased for each element by adding a real constant (3.6327). 
It is not allowed to be greater than 9.999; this is done by subtracting a 
real constant (8.4153) each time § becomes greater than 9.999. 

Next an exponent x = (1/10+J/254+0.62832) is formed. Next the num- 
ber f « ot is formed. Next the number | (fractional part of J) is 
obtained. Finally the relative error is determined as yx, where r is 


the preselected maximum relative error. Then 


rie ~ ee x Yr X ay 

The algebraic sign is given by the first digit in the decimal ex- 
pression of ] . If this is even, the sign is positive; if this is odd, 
the sign is negative. 

A special treatment is used if an = 0, for which ai =a 2 , the 
sign being determined as before. This has obvious disadvantages if the 
average value of the elements in A is small, and it is suggested that 
future work make an appropriate modification to the present procedure. 

In retrospect it appears that perhaps a more truly random variation 
could have been attained by using 2 = (fractional part of 100 JS a 


say. Furthermore, it would have been desirable to incorporate such ran- 


domizing influences as the date, or local time, or job sequence which 
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should become available internally to the computer, but which were not 
available at the time this work was done. 

For convenience in comparing values of C* with C, to see how good 
the recovered matrix C* was, most of the compliance matrices C used were 
diagonal matrices. However, there is no loss in generality as may be 
seen from the following argument. 

Starting with the matrix equation CU = V, with diagonal C, and pre- 


multiplying this equation by an orthogonal matrix M, it can be written as 


mcm'Mu = MV (19) 
where the insertion of vey in equation (19) between C and U does not al- 


ter the equality, due to the fact that ven = J. 


Now let 
mcm! = @ 
MU = U 
MV =V 
equation (19) can be written now as 
cU =V (20) 


Equation (20) has the same form as original matrix equation CU = V, 
but with the difference that starting with C (a diagonal matrix), the 
new symmetrical matrix C is no longer a diagonal matrix, and instead all 
its elements are generally non-zero. 

Thus a general problem CU = V could be written in the particular 
form CU = V, with diagonal C, if only the orthogonal matrix M, which dia- 
gonalizes C were known. C = MGM, U = MU, and V = MV. 

The particular validity of this argument was verified by two types 


of numerical experiment. In the first type, nondiagonal C's were used, 


errors were introduced in U and V = CU, and the recovered matrix C* was 
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successfully compared with the original C. In the second type, an ortho- 
gonal matrix M was generated, and matrices C, U, and V were formed using 
M, diagonal C, U, and V. Then random errors were introduced in both sets 
of data (U, V and U,V) and C* and C* were obtained. Comparison of Mica 
with C was quite satisfactory. 

From theory the minimum number of independent pairs of vectors F and 
D needed to solve the problem is 4. This is because in the compliance 
matrix C there are 21 unknowns, and 21 independent equations must be form- 
ed. Thus would require four pairs of 6-vectors F and D since the notion 
of 3% pairs is excluded. However, the particular methods employed in the 
programs seem to demand a larger number. The minimum number of pairs of 
vectors needed to obtain "good" results was six. 

In one case, a problem with 24 pairs of vectors was solved twenty 
times, using a maximum random error of ten percent. The errors were in- 
troduced in a different manner for each of the solutions. From the twenty 
compliance matrices obtained in this way, an average compliance matrix 
was obtained. This matrix looks ''very good", when compared element by 
element with the original compliance matrix. This indicates that if the 
value of n is so large as to overload computer memory or as to prove in- 
convenient for loading data, the U and V matrices may be conveniently 
partitioned, a number of “smaller problems" dealt with, and the answers 
averaged. One could regard the case described as forming a single pro- 
blem with n = 480 which was actually treated as the average of the re- 


sults of twenty problems, each having n = 24. 


37 


SECTION 8 
CONCLUSIONS AND RECOMMENDATIONS 

There does not seem to be any direct way of testing the validity 
of the theory and the computer programs which implement it other than by 
taking actual examples. These examples were constructed by assuming 
values for matrices C and U and generating matrix V = CU. Then, in most 
cases, errors were applied so as to get modified matrices U* and V*. 
Then, from these matrices, considered as input data to the program, the 
various procedures were used to recover an optimal value, C*¥, which is 
then compared with the original C. No useful or convenient method has 
been found for assessing the accuracy of the recovery other than by com- 
paring term by term, and, to simplify this task, in most cases C was 
assumed to be a diagonal matrix. It was shown in Section 7 that there is 
no loss in generality in doing this. 

One of the things that might be expected to influence the ease and 
accuracy of the recovery is the relation between the eigenvectors of the 
matrix C. Tests were conducted, at the one extreme, with diagonal C 
having identical diagonal elements, and at the other extreme, with dia- 
gonal matrices having diagonal elements differing by a factor of 3125. 
In all, fairly extensive numerical work was investigated and various 
maximum percentage errors were introduced. The output produced as a 
result of these studies is much too volumnious to be reported in detail 
and it has not been found possible to reduce the volume of information 
in a satisfactory quantitative way. Accordingly, it is necessary to 
restrict remarks to those of a general nature which will be made in the 
following. In general, it may be concluded that the theory and the pro- 
gram have been proved to be accurate and effective, and in combination 


they provide a useful tool for analysing the data from certain types of 
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structural and mechanical experiments. 


The following matrices have been used as compliance matrices C. 


from which the matrices C* were recovered. 
9. O 0 0 0 1. O 0 0 0 0 
0 Sls 0 0 0 0 2 0 0 0 
0 0 : 0 0 0 0 Se 0 0 
(1) ; (2) : 
0 0 O° 12 0 0 0 0 0 4. 0 
0 0 0 8. O 0 0 0 0 5. O 
0 0 0 Que27 | 0 0 0 0 0 6, 
1. O 0 0 0 0 fe O 0 0 0 0 
uy FO. 0 0 0 0 7. O 0 0 0 
0 0 100. 0 0 0 : 
(3) (4) 0 7 0 0 0 
0 0 0 100. O 0 0 0 0 7. O 0 
6 oO 0 0 10. 0 | 0 0 “Om Cee 7a 0 
0 0 0 0 0 ly | 0 0 0 0 0 ae 
> 
5. O l. O 0 Le 1. O 0 0 0 i 
0 5. O 0 2. 0 2° 0 0 Leavo 
, : ings 0 : : 
(5) 1 0 5 0 0 (6) 0 3 1 0 0 
0 0 0 Seee0 0 0 0 le Aa ele ae 
0 a 0 5. 90 0 i ies 0, i 25. ee 
pss Deal 0 0 De i oO 0 0 lL. Gabe 
lee 0 0 0 
0 Oe 0 0 
0 0 ho 0 0 
(7) 
0 0 OO - 125, 0 0 
0 0 0 0 625. 0 
0 0 0 0 085125. 
OUTPUT No. 1l Gives the desired compliance matrix that satisfies 


criterion one. With this procedure the types of 
matrices numbered 1, 4, 5, and 6, have been recovered 


quite successfully using up to 35% maximum random 
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OUTPUT No. 


OUTPUT Nos. 3A,3B, 


OUTPUT No. 


2 


3C 


3D 


errors. Matrices numbered 3 and 7 were recovered with 
up to 54 maximum random errors, and matrix 2 was 
recovered with errors up to 25%. (Where no errors 

were introduced, each original was accurately recover- 
ed). 

Gives the desired compliance matrix that satisfies 
criterion two. All the matrices were recovered same 
as the originals when no errors were involved. Matri- 
ces numbered l, 2, 3, 4, and 5 have been recovered 
quite successfully using up to 35% random errors, 

and matrices numbered 6 and 7 with up to 25% random 
errors. 

Gives the desired compliance matrix that satisfies 
criterion three, using the procedure that obtains the 
Square root of a non-symmetric matrix. All the matri- 
ces were recovered same as the originals when no errors 
were involved. With matrices numbered 1 to 6 this 
procedure has been worked successfully with up to 35% 
random errors, and with matrix numbered 7 with up to 
20% random error. 

Gives the desired compliance matrix that satisfies 
criterion three, using the procedure that takes square 
roots of two different symmetric matrices. Matrix 
numbered 1 is the only one that has not been recovered 
same as the original when no errors were involved. 

In the light of other results obtained, this failure 
to recover matrix 1 to a very high degree of accuracy 


is unexpected and it has not been satisfactorily 
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explained. Matrices numbered 1 to 6 using up to 35% 
random error, have been recovered quite successfully, 
and matrix numbered 7 with up to 25%. 

OUTPUT No. 4 Gives the desired compliance matrix that satisfies 
criterion four. Using any value for BETA, when no 
errors were involved, all the matrices were recover- 
ed same as the originals. When the error was increased, 
better results were obtained for small BETA‘'s. The 
same remarks given for OUTPUT No. 1, applies to this 
output. 

To have a better idea how these results look, see Appendix C, where 
the computer solutions for a particular problem with zero and ten percent 
random errors, are given. 

The work presented here is just the beginning and not the last word 
in this new area; more work and investigation should be done. Some sug- 
gestions are given below: 

1) Solution of equations that satisfy criteria 1, 2, and 4, 
can be investigated using the procedures given by Bickley and McNamee 
[4], or possible modifications to take advantage of the fact that matrices 
A and B are equal. 

2) Investigate how large maximum random error percentages one 
can give to matrices U and V. Give to these matrices different maximum 
percentages of errors, e.g., to matrix U, 5%, and to matrix V, 15%. 

3) There exists a fifth criterion given by Brock [1], that has 
not been investigated yet. 

4) Devise a new iterative procedure for the equation that 


satisfies criterion four. 
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5) Devise a convenient criterion which measures the "overall" 
accuracy of recovery of a matrix C. At present, one merely compares 
element by element. 

6) Employ advanced and sophisticated statistical theory for the 
analysis of errors and the accuracy of recovery. This would seem to 
imply constructing an adequate set of definition and the use of probabil- 
ity spaces, etc. The present analysis, in which general conclusions are 
attempted from specific cases, leaves much to be desired. 

As things stand now, the engineering analyst has, in this theory 
and program implementation, a four edged tool which can be used for the 
analysis of experimental data. The decision of which criterion yields 
the "best'' result is one which must be made on the basis of experience 


with the physical problem. 
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APPENDIX A 
INSTRUCTIONS FOR USE AND LISTING OF PROGRAMS 
A-1 Preparation of Data Deck and Final Assembly of Programs. 

The preparation of the data deck is the same for programs CECI l 
and CECI 2, but has some alterations depending on whether MODE = 1 or 
MODE = 2. 

In what follows the instructions of how to prepare the data deck is 
given for both types of modes. 

The first card is bbbbbbb6NNbbblbb....., where b = blank, the 6 is 
the number of rows, NN is the number of columns, and 1 means first matrix 
to be read. When MODE = 1, NN = 06, is the number of columns in matrix 
C3 when MODE = 2, NN varies from 06 to 50, and is the number of columns 
in matrix U. (Thus U is limited to a maximum of 50 columns). 

Following this card, comes the elements of matrix C or U as the 
case may be, punched in five entries to a card, in FORMAT (5F14.6). 
Entries for the first row are thus entered until all elements for the row 
have been punched. Then, starting with a new card, entries for the sec- 
ond row are punched, and so on, until the entries for the sixth row have 
been entered. Each row starts with a new card. 

Card bbbbbbb6NNbbb2bb....., follows. Then the entries for the 
matrix U (if MODE = 1) or matrix V(if MODE = 2) are punched in the same 
manner as before. 

Example: Punch cards as in Figure A-l, for the following 6x6 diag- 


onal matrix as C matrix. 


e 


COCOONS 
OO OWO°O 
ooFroo°o 
OmMod oO Oo 


oOoOoOO0O OF 
oko koekoeke) 
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bbb6.0bb 





cementite 


Sixth row B BB B bbb5.0bb 


—_o 3} 
Fifth * 


Fourth row —~S\, 


Third row 








B BB bbb4.0bb 







Second row B B bbb3.0bb 


First row 






First card B bbb2.0bb 






bbb1.0bb 
bbbbbbb606bbb 1bb 





NOTE: b = blank, B = 14b's. 


FIGURE A-1 


Note that element C(1,1) is punched between fields 1 to 14, element 
C(2,2) between fields 15 to 28, and so on. Note too that each row has 
two cards because there are six elements for a row, and in each card 
only five consecutive elements can be punched, 

If more than one set of data is used, the second set comes after 


the last card of the first set, starting again with the first card 


bbbbbbbNNbbblbb....... 
and so on. 
In order to stop the program, the final card of the data deck is 
punched as follows bbbbbbbOb6bbb6bb....... 
The complete program deck includes: control cards, main program, 


subroutine package, and data deck. The assembly of these components is 
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shown in Figure A-2, where the first control card is from the CECI l 


program. To use program CECI 2, change the word CECI 1 to CECI 2. 






Data Deck 


//GO.SYSIN, ,DD,* 






















SUBROUTINE PACKAGE 





cc 




























N 


AA 


//FORT.SYSIN, .DD, * 


//CECI1, JOB,,'9725002', 'VILLARAN', MSGLEVEL=1 







a 
la yi 


CARDS = 





NOTE:, MEANS BLANK SPACE 


FIGURE A-2 
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Am2 Listing of Program Cw! 1. 


The cards of prozram CuCI 1 are numbered from 0000 to 0159. 


CORA FIN OM SOHO ANMOA PIA OH OROQGANMFINDE GHOMAMFINOMH DHOKHKNUMH TN Of 
DO DODD DODO O Ate eta ed et eR OUNUN ANIA NN MOM MMO MMAMM MMS SSeS 
DO ODODODONOCDOCOCDCOCOCOCDDVOCOCOCOOCOOCODOWDOOOCOOCOODOCCOCOCH 
DODDODMVOCODOOCDOCOCMTOCOOCCOCDDOCVOOCVOOCOMOOVOOCOMODOMODOOOCOOO 


= WW Ld - 
o WJ WY i | CO ner & e 
~ > aD iwi 2 © a ry) . Y o>d< Vp) 
N me ESL {uJ o 9 D> an) Om ete em | LL 
Cs COOl OT HNIONOEAERK ON KPWNHZKELD M+ ac >< mm) 
om > Ae Se a NWOU TdF CeKe ewe! ae bm ee LES OS tee 
_ << PS Wk eam IY GMb UW © mod OT OT fe 
0 << TO wd KILO OF me Fr ome Zit 0 Fer ate 
oo Met OmSus ee ae Ww eWNOedOQ irae & Lf 
= Ym RAeis i Terete#eZe OO Owe = 
LAA] ZS>5 OF =a NeceagrOWw llaee Yea OT OO 
ww AJ Ww fa dys Seam ea ORT ~ ¢ z—_— YM ANOW <{(Abe LLL oF 
(a he COrO ww em eo =—_ Ww wee 2 Gow 
hk oJ AZO AMweae MUN eH OOM df iti bh eT ke 
1 a) te YK SF DOUSKSZEOCINOZO -& ~ee Se 
—O MW TWYOwWwCA e OCALA rRaAaeaAaA Ow Ta OwWSCeY 
© DS ewWs C2 ST ehedus <= CL LL! Co CMOaere~wWw WO 
we oc eft et NYHOWLEL ee De * Ce Wet SelM 
om FO ener uvOQO OO wf mDACT et EAU mnt 
aaa WES Ieee A e¢ TOW Aa C2> <I LL > p= bY) 
— © S2MIWAMRMNOwOdWwWreE dT OZT eZ2rwoo miedo << 
o> wu — STW F< Whe OO me Om a OG: CCl & 
om ef WX TYNDWeE FE OLD KW MDM Cwikft mt em OY OO om Li 
QO om mee — OD PW ede Ss erm Ind etOD wat Co 
er an oh | 8 eg ZzZMWSS>O0O AZ2aocaw JO Ware WW Wl 
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APPENDIX B 
SUBROUTINE PROGRAMS 
B-l Description of Subroutines. 

SUBROUTINE ADDEM (A,B,C,D,M,L) provides the result of A = B+DC, 
where D is a scalar constant, and A, B, and C are double precision matri- 
ces dimensioned (M,L), where M and L are not less than 3. With D = 1.0, 
A = BIC; with D= -1.0, A=B -C. 

The matrices B and C from the calling program must have identical 
order, otherwise a message is printed and the execution is halted. 

SUBROUTINE AINVM (A,B,C,D,L,N) forms the inverse A= 5 of input 
matrix B. Matrices A, B, C, and D are double precision, with A and B 
dimensioned (L,N). Matrices C and D, which are used as scratch matrices, 
have dimension (N,N), where L = M+1, and N is not less than 3. 

A message is printed if the input matrix B is non-square and the 
execution is halted. 

SUBROUTINE BIGCM (COL,Q) forms the 21-element vector COL from the 
6x6 matric Q, according to the scheme described in Section 3; if the 
input matrix Q is non-square a message is printed and the execution is 
halted. 

SUBROUTINE BIGSM (BIG,P) forms the 21 x 21 matrix BIG, from the 6x6 
matrix P according to the scheme described in Section 3; if the input 
matrix P is non-square a message is printed and the execution is halted. 

SUBROUTINE BLANCA (CO,BE,N) forms matrix CO from matrix BE by intro- 
ducing random errors in the latter according to the method described in 
Section 7. The arrays CO and BE are dimensioned (7,N). The maximum 
relative error, denoted by the symbol r in the exposition in Section 7, 
is governed by the computer variable PERC, which is one tenth the relative 


error. If PERC = 1.D-2 = .0100..., then r = .100..., and a maximum of 


2. 


10% error is introduced. PERC may be changed by changing card number 
0313 in the subroutine. 

SUBROUTINE JACVAT (B,N,NOYES, RR,SS,EIRU,NDIM) gives the eigenvalues 
and eigenvectors of the symmetric matrix B. If B is non-symmetrical the 
subroutine first makes B symmetric by forming p+pl RR is the spectral 
matrix that contains the eigenvalues of B. SS is the modal matrix of B 
that contains the eigenvectors. EIRU is the square root of the elements 
in the diagonal matrix RR. N is the dimension of B. NOYES is 0 or 1; 
if 0, no eigenvalues will be recovered; if 1, the eigenvalues will be 
recovered in matrix RR. In the application made in this thesis, NDIM is 
the same as N. All the matrices and variables are in double precision. 

SUBROUTINE MAUSS3 (N,EP,A,X,KER) is a double precision FORTRAN IV 
language subroutine which inverts a symmetric matrix. A is the matrix 
to be inverted, X is its inverse. Both have dimension of (N,N). EP is 
a tolerable error, taken as 1.D-20 in the present work. KER is an output 
flag having the value 1, if all has gone well, or 2, if there has been 
any difficulty of execution. 

SUBROUTINE OUTPU (R,B,U,V,T1,T2,T3,T5,LL,L) has inputs R, B, U, V, 
where R is the desired compliance matrix . B is the desired stiffness 
matrix, and U, V are the data for the problem (such that RU = V, approxi- 
mately). It prints the matrix R, and computes and prints the three 
criterion numbers OQ; (i = 1,2,3). Matrices T1-T5 that appear in the 
calling statement are scratch matrices. LL and L are dummy variables 
used in the dimension of the matrices and are related to the number n of 
vector pairs used. 

SUBROUTINE PRINM (A,L) causes the (6,L) matrix A to be printed, 
where L is not less than 3. The first line printed contains the fol- 


lowing information: MATRIX NUMBER XX, YY ROWS, ZZ COLUMNS. This is 
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followed by a print of the matrix, line by line, in FORMAT (1X, 8D12.4). 

SUBROUTINE RECIP (A,B) A and B are diagonal matrices (6,6) and A = 
pial where B is the input matrix. 

SUBROUTINE REDEM (A,IDENT) reads into the computer memory from the 
data deck which is punched as described in Appendix A. A is the "name" 
of the matrix and IDENT is an identifying integer to assure that the 
card-group which is read is indeed the one desired. 

SUBROUTINE SQRERT (A,B,C) solves the matrix equation that satisfies 
criterion 3, using the procedure that twice takes the square root of a 
symmetric matrix. B = uur, Cc = wi are the inputs, and A is the desired 
compliance matrix. 

SUBROUTINE SQROOT (A,CC,G) takes the square root of a non-symmetric 
matrix by iteration procedure. Matrices CC and G are inputs. CC is the 
approximate square root; G is the matrix whose square root is desired, 
and A is the desired square root. 

SUBROUTINE SQUAM (A,B) forms a 6 by 6 matrix A from the 21-element 
column vector B, arranging the results properly as implied in Section 3. 

SUBROUTINE TRACM (T,A) forms the trace T of the 6 by 6 matrix A. 

SUBROUTINE TRANM (A,B,LL,L) forms the transpose of matrix B(7,L), 
obtaining matrix A(LL,6), where L is no less than 3, and LL=L+ 1. 

SUBROUTINE ULTIMA (A,B,C,MM,M,L,LL) forms the matrix product A = BC. 
A(MM,M) = B(MM,L) x C(LL,M), where MM, M,L,LL, can be any numbers not less 
than 3 for which the matrices are conformable for multiplication. A 
message is printed if the matrices B and C are not conformable and execu- 
tion is halted. 

SUBROUTINE ULTMB (A,B,C,LL,L) multiplies the square matrix B by a 


column matrix C, to obtain the column matrix A. A(LL) = B(LL,L) x C(LL), 
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where L is any number not less than 3, and LL=L+t+1. 

SUBROUTINE OUTPUT (R,B,U,V,T1,T2,T3,T5,LL,L,BETA) is the same as 
SUBROUTINE OUTPU. The only differences are the inclusion of parameter 
BETA in the calling statement, and the printing of four criterion numbers 


instead of three, o%&;, (i = 1,2,3,4). 


B-2 Subroutine Listing 


CARDS NUMBERED 


SUBROUTINE FROM TO 

ADDEM 0160 0182 
AINVM 0183 0208 
BIGCM 0209 0233 
BIGSM 0234 0288 
BLANCA 0289 0336 
JACVAT 0337 0481 
MAUSS3 0482 0532 
OUT PU 0533 0553 
PRINM 0560 0576 
RECIP 0577 0594 
REDEM 0595 0616 
SQRERT 0617 0638 
SQROOT 0639 0692 
SQUAM 0693 0727 
TRACM 0728 0744 
TRANM 0745 0763 
ULTMA 0764 0791 
ULTMB OS) 0814 
OUTPUT 20000 20360 
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APPENDIX C. 


SAMPLE PROBLEM SOLVED BY COMPUTER 


C-1 Sample Problem Using Program CECI 1 with 0% and 10% Maximum Random 
Errors. 


The nomenclature used in the output of the program is as follows: 


il 


MODE 
Matrix 
Matrix 
Matrix 
Matrix 
Matrix 


Output 


Output 


Output 


Output 


Output 


Output 


1 

P 

BE 

EB 

U 

V 

No. l 
No. 2 
No. 3A 
No. 3B 
No. 3C 
No. 3D 


Data will be generated in the program. 
Original matrix C as read from cards. 
Matrix U as read from cards. 

Matrix V from V = CU. 

Matrix BE with random errors introduced. 
Matrix EB with random errors introduced. 


The desired compliance matrix which satisfies 
criterion one. 


The desired compliance matrix which satisfies 
criterion two. 


The desired compliance matrix which satisfies 
criterion three, usingthe iterative procedure for 
Poa the square root of the non-symmetric matrix 
vvtuu~. 


The desired compliance matrix which satisfies 
criterion three, using the iterative procedure of 
meni the square root of the non-symmetric matrix 
UU ‘ 


The desired compliance matrix which satisfies 
criterion three. It is the average of outputs 
3A and 3B. 


The desired compliance matrix which satisfies 
criterion three, using the procedure that twice 
takes square roots of symmetric matrices. 


The printed output for the problem with zero error appears on pages 


77 to 


pears on pages 


82 inclusive; that for the problem with 10% maximum error ap- 


to 88 inclusive. 


76 


MATRIX P FOLLOWS 


6 COLUMNS 


6 ROWS, 


ly 


MATRIX NUMBER 


7O00D 02 


Oooo0o Ww 


OOCOOO 


000D Ol 


OO00C0O 


OO000°O 


200D 02 


OOO HOO 


COooo0°0 


OCOD Ol 


OOMoo eo 


Oo0000 © 


MATRIX BE FOLLOWS 


29 6 ROWS,24 COLUMNS 


MATRIX NUMBER 


CIOIAIN eOUOIOUINN AION ION 
ealelolelalelelal dalolajlalelelalea 


COBoaaooanoonoonaoaeaoado 
OOO OCOMVOCOTC OOCOOO0OOC00°O 
NSOHOHOOMOODOOOOOW 
CNM NINGAE NO DO AIAN MINS OID 
WO FUNG M OP COM DH Not FN ot 


de: ay eas as aes camer aaa 
| i ! 


OVODOCOO0OO 


AOCCCOCO006 
OOODOCOOOO00°O 
COrDBOW ODOODOOHOON 
ODOORS TNA ANAM HMO AA CO 
MAMAN A MOaeEMMOMODAHNDOONM 
ecoereeeoeeeoe eee e ee o 6 


thr ir asiiasiainaslis acide’ wim 
| 


ONNANAGA! 
Oo000° 


(anny nl ee Tay 20 ae) 
OOO0C0OO 


Nat QUuaQian 
OOOoO°O 


QAQ0aG 
OO00oO 


MICU QUAI QUAN QU ONAL 
alelelelalalelelejlaiele) 


AQAOOOQAQGHOAOGSO 
COOoo0d0d00o00n0o0co 
MOWMWMOneOOOOOO DOO00®D 
TNMOPEPMMNOCEROM WOME 
OS ed ON ed OP £0) O SPUD COLA © ot 0 TM od CH 
@®ee?eee#eee@eeeee#e¢@e@ee##e@ @ 
alain acai ac inl i? 
| 


AIAN QUINT O QIAQUIQU QUAI A et QU AIA CY 
elelealelelalolealelalealejlelelaleal ol a 


OOOO OOO08S8O0O0G00000 
OODOCO ODOT OOOO OOCO000O 
OO DOOOCOOCO HOOMmooodU0oem 
BRAT HONS DINO a O et MLA 0 
OP mt OP LA OND COCK AI TLL SN 
ecoeocecueoeceeaeoeeeoe eee eo 
5 ee ee aa 
! 


NINE AINNNNA NOAA et OU 
OOOO CO0O0OO BOQOO000O 


QAGAOOGOGOGRAaA ASOO00006 
OOOOOO9000O OBOO0O000O 
OOOooOMmoo0oo0 OVCOCOMNMOO 
DMAP AAR OU ONS UO O 
SL PD UD et Diet O CULV EOLA OO) SHO 
eeoeeveeceeeeoeeeceee een 
ee a | 
! | 


QIEOIAIAIOIOI QUAI AIOIEAIOAI OU A teu 
OOOO V9OO9OCOOO00000 


OOOOODO0000000000 
OOODOCOCOCOOOO0O0000 
OOOMOOVWCOOOO00 000 
AQUINAS SH aNa AANA OS 
md POD OOF) ed 0) rd OO LT NU SO 
eocoeeeveeevees eevee 6 o 8 
ae ek ee 
' 7 ! | tt 


NNANNNNANAN Ne ONONAIA 
OOO OOVO VO VO O09 V0900 0 


ADD OONOGSO0O000000000 
OOOCCOCOOCOOOOCO0O0O0000 
OLTOOPFROOQVOCOCO9O0O000°O 
NOANOIAEIN SF HOO RON 
POW SNMANMA PMO SNM AO MAA 
eeeoeoeveoer eevee eeeeee @ 
i iid (nai WaalmMMO ini 


NNO et AIN NI AIAIAIOUIAIA Ht QuQluaia 
OOO OOOO O 9000900000 


OAGBOQOO0O000000000080 
elalejleolelolelelelelele)lelelejelela) 
OQVOPrF ODOC OCOCONMOO 
MmNOODT NEAT TOMOTOr 
DOA OD Pin Oat Om NA PUNO SP 
eoeeveeoveoeeeoe eee 0 6 © @ 


+ tats iniiiad Gali ba Dion’ ~via Cm 


900D C2 


77 


oo 


MATRIX EB FCLLOWS 


6 ROWS,24 COLUMNS 


3¢ 


MATRIX NUMBER 


COMM MAILE MMMM MMA LMM 
OCVOOTOOOCOOCOGTOO080 


Ooeaghooonckaoanecoaaceac 
DOTMDOTDOOMOCACCSFTOOW 
Set LAE OPM NOM ONC FOOD ff O 
WWM DOMAST MNOM MOM HOOMINAN 
SUF OPM HOOT ON SHAWN 
®eeeeee#8e#teee8eee8s@¢8%®©6.e6%©6h6mc8 Ch 
Og a OA apa al 
l | { 


MAA AAO ANAM MOMMOLOS 
COOCKO TOOCOCTWOOCOSO 


OCOOQaAa0 COOCOOCO00G 
OOMTOO GCOCOOOODMOr 
ONO LTO OOONOASTOLTOMM 
MNMM-ON NMNODOMOODOSLDO 
DO AID ENON OC et et LTC et 
eeeoeeeeeeeeeeetke#eee ee 
COO OS OS a aaa 
i | | | 


MAM OF FO MAIC CO 
ODOOOCOOOoOCS 


OeCOoocoaecoeeadhsd 
mONMNDONDOCVO0CO 
DDAMROMWMNOODANO OCOOMMA 
HONKONMOMOMAOO NOTOO 
CN tN ed ed etd OO OOD SP OOO OD Sf mt 
eeeeeee#8e#@6¢e¢eeeee28e0e80e?08 @ @ 
Se oan 
l i 


AOO PHL HAMA POONA S St 
DCOOVCOVCOOCOCOVCO0O00O 


Ooo ooeon fSooaeeaoeeasS 
OOTONODCONOOOCOWOMMO 
“N= FNOOOONODOOTOUM 
OMmOAMm SFOAINAOP NODNANANW 
0 erred ered omed OPS OY ered SS md ION SE OD pt ed oe 
eoeeeeeee%eeeee?e 8 8@® @ @ @ 
‘ it : [ 


ONOOS 
Ooecdco 


Oooo es 
OOON0e 


AOANUSOMONMON NAMOMOMOM SOS 
OOVOOV OOOO WBWOCO0O00O 
OCOooooaoaoaoo aooaaeaace 
ODOVUOODCCONO WDOCOODON 
NAMOSTOCOWO NOWUDONONW 
TMDMMAINOWM MWANMODONAO 


FOAM ANEH ONODNSENRAIN 
eeteeeeteeeernretrteeeve @ 


| | 


MOF LOAM MOMMA OOM MAMO MAS 
OVO OVOO ODO GVOOO000O 


QDAOOGOAOOOOODOO0O000 
TOOK OCOVOWOVOVOVOONO®D 
TOK OOMONTNTOOWA 
TOP HEH NNDKNOWMOS AO 
OO CF) red red nd UF) ceed ed PO UD ed WO eed ed 
®oeeeeeeekeeeeeeeee?ee ee 
OOODVVIO VO VWOOMOCVWVO0O 
' i | | i | i | 


NAMM OO COM MAM AI OO) a AICO ON 
OOOO DOVOOVOVOVOO9O0O 


ADOODODOODOO0OO000000O0 

ODOSONNODODOOOOOOOCO9O0 

DINSOSFTOWMMNODNODNOLEH! OH 

ND OHOS OM HOODS Sete ovdned 

Om ST OST OMBMAINSIENT OOO 

oeeeeeerekeeee7#e#2»eeeesee%tk eee e 
| [ ! itl 


NOD HA LOAN AO OOM OAIOS $F + 
OOVDVDOVD VDOC DWDVOMDOOCOOC CO 


7100 C3 


OOdOOOOOOOOOOOO0O0000 
OOV0V0OrF OOOO 00000u00rm 
MDBONM FOMNDWONQOOMOANON 
Mm ONO DAT OVONOVOMNADGM 
CD ONY cd dS ONS C0) reed Pe od OA) eed (O) raed eed od 
e*eeeeesrtsee?¢??eee0ee0e0080 0e@ @ 
DOODVOOCOVVODVOVOVVVMO00°O 
! | 1 | f l 


C.0 


ERROR IN BLANCA 
ERROR IN BLANCA 


~~“ 
C2 


0.0 


MATRIX U FOLLOWS 


6 ROwWS,24 COLUMNS 


4+ 


NUMBER 


MATRIX 


QNINN RON ONINIONIONN Net 
OCOCOSCOVOCOCOOOOCCOC 


CoOOCOoOCOOOCOC OAS eeaw 
OOWVCTOCODOOVOCVOCOCWOUW 
NPOHKHOMOOMOOCCOMoOoOor 
ONAN AE NN ODHINMIN SONS 
WO FUND DWE OP MP CON at FA ot 
eeeeesrseeeee%s?8eeeeeeeee 
in ae ee! 

| 


ONAANANAN 
COCOCOO 


OOAG006 
OOCOCOO 


NAN ANNAINAIAIAIAN 
OOCOSCOCCO0OG 


oOoocoocoeeeae 

OOCVOOVCO0CO°G 

OOFDOWMW DOOCCDOOOON 

ODOVOOE SF NNKOAMNM IOS OO 

—NONAKMOMORMOM ODN OM AM 

eecevweeveeeeeeevre ee ee e@ 
[ | 


NOUNANNNN NIAAA 
OCCoOOCOOCO0O0oO 


OOoooooo0ooo0o0eo 
elelelelele}l elejle]e]e}e 
NOWMOABMOOOCOOD DWOOOMH 
TNMOP-NMNODOMOW MWONIN™ 
Nat NOM MOFINMOMNOASM Hee 
eeeee#2tfe?8e#s8e?e?eteee¢e?ee e 
ONMOOOOOCOOTOOVOCOOO 

| \ Vil I 


NANAINONNAAIN Neat NaNO 
OD OCOCOCOOCOCCOCOCCOO 


OQOoeaooCcooooaoeooease 

SOOO OCOOCOCOOOOCOCOOCOCOO 

COODOOO0 OOONMOQO000% 

AUTH AONTONM HOM MOAM O 

PHAM NONE DANN ANN SiN 

eeoereeevereeer ee eee & 0 
i ' 


NIN SAINI 
OOOO OO000 


BDOoOooOoodaaooed 
ODOCCO0CC0O0O 
ODoOoOoOmMnoodod GOOCONNMNOO 
DHOOM ANAM AOA O O 
MAF DOAN HALA HH OAIAM LACH SF DP 
eeeeeeeeeeeee?eeeeee @ 


OOO OO OVO OCVOO0O00°0 
[ it | 


Ne lainy 
elelelele) 


COoOoOaoO 
OO0COoOO 


NQINAINANA ty 
OOOO0O000 


COoceceaooaoa 
OOOCO000O 


NIOAINAIONIOAIOIOAIATOI AI OAT A at ON 
elelelel ee] eiealeloleleal@eleleleal® 


QOOOOOO0O000000000 
OOOO VWDOOO9VO9O9O000O 
OOONMODOOGCOOCVOO0O0O 
SOWA ST SAHIN DONNO 
OPN DS OOO et (ot FOO ILA ONIN SO 
eeeeereeeereeeev ee eo @ 
Be oe Oe 
' I | [ ! | 


NNANNNANN ANNAN CUNONAIAN 
ODDO ODVOCVOO9OODOVO000O 


BOOODODOODOOOHOA0Q0O000O 
ODDOCVWDVCVOCOCOCOCOCOCCOCOO 
OTVONARODDOVVOOV0O0O 
NDANADOINMINA ST HODDONAMD 
ROMFPNUOANMALOASNNMO Ma 
eeteeee?#eee#ee#ee@eeeee e@ 
Bo 8 O00 ann 
1 


NNO MAIO OIA AION ON ON 
ODOT CWO OOOV OOOO OO000 


QDOAQDOOOOOODOOQADAONOO0OO 
DOVUDVCOOVOVODWVVOOVOVOO 
OOVOFrFVODCDOCDCDOCOCONMOO 
MmNOODINADTSTOMOTOFM® 
MOANA AWW SU OAD FINO SS 
@eeeee?@eeeeeeeeeeeeee e@ 
DOOVOVVDOVOOCCVOVOCO00 0 
| ! | [| | | 


9000 C2 


MATRIX V FOLLOWS 


6 ROwWS,24 COLUMNS 


-~O0.17100 C3 


D9 


MATRIX NUMBER 


AMM MONS MOCAM MAM MAIN ST OAM 
GOOOTOCCOOC COCOOCSO 


SCcanooocaeoanooooaooccona 
TDOOCODOOMCCACOOtT OCW 
aL OMA Ost ONO FOO SS O 
WOOD SF MOOQMOM HOODMAN 
FAP St SoM TODS OA es HtiN 
@®*eee?ee?eeoe¢e¢e?0e6e¢¢0e¢e0e0¢06h¢60668e@0 «@ 


OOOOOCSCOCCCOOOO0O0°0 
j j j | 4 
KN OA CAH CO 
Cacnoondc 


caceoc 
COMFOO 


NAIC MO CVC EO HO SF ON 
COCoCcogoco0o0o 


CSAGCOGHOQ000QG 
OOOCOVOODMO 
ONOFOCH COOCONA S£OLSonm 
OMMM-ON MDOOMADOVOSTDO 
DAE CACAO © PB tt Pe Pm et HLA Hs 
e*eoeetees?#?t¢t?eef © @ 6 8 OO 
OO OS aan? 
| | { 


AOA AMF FOAM NMMOM 
OOOceococqao0co 


AMOQOGOOO000O 
M-ONnODONOOCOCOO 
DOH4ROMMWOONNG OdIAinn 
AOMONMMOMAMAO NOPTOO 
CI ot A ed et eet mt OO POO HOO SS 
®eee?eeoeee?#ee0e00 0 ¢¢ 406 80 @ 
SO O0C OOO TT 
j 


MOOS PTHOMOAMGSCOMONNT SS 
OOCOOCCOOOC OC O9OO0O0C Oo 


OooooceceogcaaacoooooQ 
OOTONCOOHOCOCOCOOHOAO 
HN FNUCOOONDDOOLTOWM 
MOM PFORNA OP NODNNAIWY 
CO red ad weet OS OO) eet a met ONIN SF OD ad ed 
@®eeeseree%rv#H#eee2502#0e?#28e#%@e?¢8 e888 @ @ 
O00 08° aaa a 
j 13 | 


MAN SOM AICO 
OOOOCO0O00 


SeaOoeocooeg 
OODOCBCOOO00O 
NAMOTDROWO NOWVDODNONW 
TWVTBMWOAAWOW WOWO ONO 
TAS AICI AAI ONDA NINA 
eeeteeqe%e#%see%@¢0¢¢e0¢0¢ 84 ¢@ ¢ 


O08 Oa i. 
1 ‘ ' 


CANIMOCA 
OoOooed 


QOOaagago 
Oooo 


MMA OAM A SO SF 
OOOCCO0O00 


Hooaoaceagaoo 
OOO9OODON 


MOF PO AM AM COMO AMO 
OOCWOVO OO9OOO0O00 


OdNOO000000000000 
TOOCTOOODOOOOOONO ® 
TOK OOWON IS NPOOWN 
TOM HMM OPOMO LIMP 
CO OO) cnet peed Od LEN met weed SP OO UN ad SPO nd nd 
oeeoeveoev eee ee © © © @ @ ® 
if j jf 


EVA CA CALA FY FO FO CAPO COA AICO et OO FAO 
OOOCOCOCOOOCOCOO0C09Oo 


AQOODOOAODOOO00S8000000 
ODOTPTONNOCOOCOOOOCOCO9OO 
DWPTPOTOWOWODNOONO LOH 
OVD OA POE HOD SP PMSF OO DO 
DOM OS Om Iusysyrnynowdy 
Cr ee ee ee ee ee ee ee ee ee ee ee | 
ce meas Vee area rg 

] j j 


DASA A POM AAA OAMNA SSS 
OOODODOOOOVOOOOO0000O 


QADBDODOOOQONOANOAO00000O 
OVODOM DOTOODODOOO0OU00NaM 
MOON FOLIMWOOSANUOOMNON 
Mm OUOM OatASTODAOOWWAOM 
CONS ad mt OS NI I CO) red fm end CA) rnd C0) pred med ed 
e*eet@@e@emhmhUCcthC hhCr hh—C<C OC‘ OWCPHhlhC<“ WCHL LUC lUlUrhhUr hUlClUDF 
ae ee ee 
| ] j 


10.37 
LeMINIMIZES V-SQUARE ERROR 


ELAPSED TIME 


CUTPUT NC. 
DESIRED MATRIX IS 


L7, 6 RCwS, 6 COLUMNS 


NUMBER 


MATRIX 


TOTOCO|AN 
gmt goed 9med peed peed (> 
a OS 

AMOaaCA 
ODM OOS © 
OOOO Or 
Naw FO mA) 
eeeeee 


CCOO0oO 
i | 


WU SS et 
pwnd pond pond pend ) pend 
yu | 
Oaaaeac& 
OM ADO st 
DOANTOST 
OTOKCOO 
CAD CO) 4 CO ot 
e®e@e¢ @ @ «6@ 
0 So 
1 


OSI NSW 
oat pend ad CD) nd od 


Warren 
oad ead ond pened pnd 
| jt 

QOOAQO 

DOM OP © 
WOM st 1 
OAOM ATO 
omnd ed CO) nd OD won 
e@ee3e%ss8 @ 
ire 
1 | 


P=ILOVUN OO UN 
CO md ad ed nd peed 

ps ddl 
OOOQQ0 
OmMUhm OW 
ODWMNDO 
OM Mey DOO 
ON md COMMA 
eevee 
OOO00°O 

ee 


1.7971800-22 
225554190-25 
4.8378730-24 


lL EQUALS 
2 EQUALS 
32 EQUALS 


CRITERION NUMBER 


CRITERION NUMBER 


CRITERICN NUMBER 
ELAPSEC TIME = 


1.48 


MINIMIZES U-SQUARE ERROR 


2y 
DESIRED MATRIX IS 


CUTPUT NC. 


6 COLUMNS 


&6 ROWS, 


33, 


MATRIX AUMBER 


FOO MIWA 
pad ond omed od od 
iyrri 
CcCaogcneagds 
OWI O ™ CO 
MAOAN SOO 
OAs at a (0) 
wat mt LN OO 
e028 &@ @ «(68 
COCMOOSO 
1] 


TM TW HHA 
pod ged pmed gmed >) od 
1's | 
oe ceas 
Mmm sr NOD 
NMA DWOwO 
OretNCOC 
=f HAI OH 
eo ¢h6c8tmlC HCC 
Oood0o00co 
| | | 


LUN SO ALLA 
and pet ome CD md eed 
ae 1 | 
oceaqceSG 
fTrnrnonr 
mmNOOST 
DMWOMNO 
mt OU a OO 
e®ee?8ee?8 *e 
OQoo0odo 

| j 


tN iN ed O sf SO 
ond ome ed od et 
' | ,i) 
an Tan} an] a | eT a) 
COOnrwe 
TOONKIN 
OM OWA 
NO IAA ttn 
eee ee @ 
OOCOO00 
bf | 


LN ATLA LAL OY 
mt A et md 
: proud 
Oooo0doo 
Oo re lh 
raoorm os 
ODOM OQeHN 
Ph tO OS et 
eee @ @ @ 
OO Oem 
{ 


ALAA ST 
© mt et pd od ot 

ripe 
Oogoano0ona 
OOoOnwTMm © 
Ost rr Al 
OOoOedd0OM 
OD pe CO ot ot ot 
eo e@?e 6¢@ 0 ® 
ST al@elelel® 

| | 


1.006397D-22 
1.975143D0-25 
€.67537T710-24 


1 EQUALS 
2 EQUALS 
3 EQUALS 


CRITERICN NUMBER 


CRITERICN NUMBER 
CRITERICN NUMBER 


0.96 
MINIMIZES TRIANGLE ERROR 


DESIRED MATRIX IS 


ELAPSED TIME 


OUTPU NO.3A, 


80 


6 ROWS, 6 COLUMNS 


529 


NUMBER 


MATRIX 


Sot Te 
pmed ped od pod =! (>) 
Pye 

eooeceas 
ODMDAKnaAO 
ett 0 © 
eknnkamkomhe) bed 
aA OA! 
e@®ee#ee® 

| 


reese Pete) 
ot ped od od CS) ed 
rtd | 
ca ea | ba en of ey fet 
WC OAIO 
OtaHror 
DOIN FOW 
OND emt AID st 
ee8¢ 8 6 «8 
Oo0o0o0odo 
er | 


tA Fars + 
rat od ed (om od 


INAS LALA 
Od ed ed od 
[ pit 
Sooges 
NOOTAWN 
NOWDOOU 
OMMOM DO = 
THs Oe 
ooo e 0 @ 
eteebecan rare 
| 


anf. 
OC) nt mt pet mt 
rit 
OOO00O 
Oecoro 
QOQOONDO 
OrmAMmO 
OD At et ti m4 
eet ®© @ @ 


OOOQ000 


2452470-22 
8.4468 380-25 
1.833622D-23 


1 EQUALS Se 


2 EQUALS 
3 EQUALS 


CRITERION NUMBER 


CRITERICN NUMBER 


CRITERICN 


NUMBER 


0.91 


ELAPSEC TIME 


MINIMIZES TRIANGLE ERROR 


CUTPUT NO.3B8, 
DESIRED MATRIX IS 


6 ROWS, 6 COLUMNS 


TC, 


MATRIX NUMBER 


WV $F POA 
ond peed med ed et CO 
prude 

QOecooqga 
NOMoyso 
DO HAA O 
00 LAT) CUR 
CO emt Ph CAs ot A) 
ee @ @ @ ee 
ODO0OO 

j ! 


TWF WN 
mad ond oed pend C) ond 
iid | { 
GOAQaAOGA 
Ph AIR 6 © ot 
MOMDMOOS 
= COO tO ~f 
AIM AI Sf OO 
e®@e@68@©6hUhrhC~CDCD 
OOoOo0oed 

if ! 


TOITNS ST 
wmnd , and med CS) ond ont 
Vid { | 
(aon) mj i ex fan] om 
M DAHO ato 
ODeowmnns 
P= Ht (AYU OM 
MAT OT MLA 
eo © @ © @ (@ 
OQO00°0O 

{ ! 


Wfratrre 
ond pond CD cand wmed ot 
' 4 (tf 
Coogee 
COMO OW O 
ONO et OU 
HOM 
COD eed ISN CF) ond ond 
e®e@¢6UcemhUOmlCU 
Oaooog 
1) 


PA ALA SO 
ome CL) med med eed et 
1 ted 
OAAQ0O000O 
OOTUUD 
TOWNOO®D 
ODE MAT 
Nd el ead ed ot 
oo ee ee 
Se a 
! 


HFSS 
OD et et ot mt me 

titted 
OQdOO00O 
OWN OM © 
OOOO 
OftOODr 
DON MUA 
eoe@e@%o8¢ ¢© @ @ 
OOO0O00O 

1 


€6-.412928D-22 
2213292710-24 
3-17C094D0-23 


1 EQUALS 
2 EQUALS 
3 EQUALS 


CRITERION NUMBER 


CRITERICN NUMBER- 


CRITERION NUMBER 
ELAPSED TIME 


OUTPUT 


0.87 
MINIMIZES TRIANGLE ERROR 


DESIRED MATRIX IS 


81 


NC .3C, 


6 ROWS, 6 COLUMNS 


68, 


MATRIX NUMBER 


WAS Po 
oand pond ond pnd ond (> 
'rd’ud 

AGBa&oa0a 


—NIFOWDO 


OWet Ono 

DON OSH 

OmaAST Orca 

e*eeee ® 

oS 
{ 


Wt fat 
pnd pond med med C) od 
Pes ' 
SQOoooo 
~ DOH DOx 
MKAtNAIO O 
ON FOO 
QUID AY mt CO 
e®*eee8 0 @ 
Gry tee 
14 


LALA AIS SP 
wmnd pond vmud CD) vad pnd 


{ irq 
AAQAAQAG 
NOT 
YO UN 4 OD 
Torneo 
US wed CS CS Pm went 
e@®e@e@¢ &@© 68 6 
Ta" 


maT TIATIO 
O mt mt ot st ot 


i'd dd 
QOOOOOO 
ODN 
ONO 
OF*AHOOV 
COUN = LANA 
eoct eee 
OoO0O000O 

' 


5-531635D-22 
1.3551570-24 
2-225371D-23 


1 EQUALS 
2 EQUALS 
3 EQUALS 


CRITERICN NUMBER 


CRITERICN NUMBER 
CRITERICN NUMBER 


MINIMIZES TRIANGLE ERROR 


DESIRED MATRIX IS 


OUTPUT NC.3D, 


6 ROWS, 6 COLUMNS 


109, 


MATRIX NUMBER 


OH OOaet A! 
OOOQOQOO 
f 

Saooee 
NORA] Aland ae 
DO HOO 
Pat OON ODO 
AFH RAN 
6° ¢©6hUC<Ci ihC<C<C HMC RHhC~«D 
OOOCOS 
\ 


OOHOxt=— 
QOOOOO00O 


| 
SaQoegge 
MM COP) OD ox 
MAOUMOS 
OOAODM AI 
ANAS es 
e©*e@e?e?¢@¢ ¢ 
OOVVDOO 
| 


=“OOaco 
OOO0O00O 


| 
[any any ae an] an 
TODMHMONA 
aA TOMI ODO 
OHM HID O 
DWN Daag oO 
e® @¢¢0 @ 


OOO000O 


OOOO 
OOO0O00 


| 

QOAQaQAAaQO 
UNO OD TON 
ONO OD 
OrMmMNOO 
an TON 

oe?ee @ 
OOd00@d 

| 


ONOOO- 
COOOO00O 


| 
OOAQQ0OQ 
PFOODO 
DONT OO 
P= COM Oe 
ad oat AIA SS 
®©ee?eee@® 


OVWO00O 
' 


—AOOOO 
OOO0OO0O0O 
| 
OdOO000O 
PTR LIAT OU 
FN CO WO ot PO 
Or ODO 
DO amt ot DH mi CF) 
eee 0 @ 


2 & 
( 


€-446506D 05 


1 EQUALS 
2 EQUALS 
3 EQUALS 


CRITERION NUMBER 


82 


&€.0641720 03 
5-876357D 04 


CRITERICN NUMBER 


CRITERICN NUMBER 
ELAPSEC TIME 


1.89 


16.58 


TOTAL ELAPSED TIME 


MATRIX P FCLLOWS 


6 COLUMNS 


6 ROWS, 


Ley 


MATRIX NUMBER 


MQOCOOO 


Soe © 


2COD C7? 


OCC a OO 


COO0OR0 


BVO Dez 


CAOOODO 
eee 8 @ @ 


OOOCO0O 


Ci 


COCE 


PT oOOuU® 


VOOOOO 


MATKIX BE FOLLOWS 


o 


van) 6 ROWS,24 COLUMNS 


NUMBER 


MATRIX 


DODD OCOCCOOCGTOO0O000 


fan [en\(=\[an\(as)am\ es (aman es any ep|[amy en} any an] any am 
CESGeOCOCOOCCOCOOCCO 
NQF CHOHCOMECOROODOOM 
ONIINNINE A ODAINMIN SOND 
IPO SINAC COMOR OE ON Ht Ad 
eoeeeee#2#?eee8e#e0@06hmhC<(aC:;2t—tCTesé~< Ow elm a le 
he SONS Stal tel apolerol=) 
| | | | 


ONANNAICNG Wet ANNNANAAAIA 
QCQOoOCeooO CoOCcoOoececo°ooo 
am | Ei = a ae ee | S| VE Ee aa am | a, Se | a a 
Oceecoqc CGeocoec Ooo eo 
OOFrDOwnNW COCDOTOCOHOON 
ODDO S NMNVNTHNONMNAMNOO AD 


SAMONHMOORMOMCOANCONM 


eeeteeoe#%ee#e8 ®#08e8 @ @ @ @ 6 @6h6UmO8hUOlU® 
OCCocooQoocooaecoqcnoceo 
i | ' 4 


NQINAIN NAINA ONION NRA 
OO@aCOoOCOCetoO ecooce 
OOO oaS eae aqaqagne 
OOOGeCICOCGeOo @OOCeeceaq 
NOWNOAONDOONG ONOAaOD 
TAMOCERNNOCACE WOW 


AI SAL = OEP ADO FINMINO HON HO 
e#eetetee@¢@#@6©8@mlUcOHmUhUCU HhUC<C OHMhC—~C MHhUC<C MMC HhC<C NHC HhC~C HhC<—~C SC! 
LO eb F ce ce X ap x en) NY «Go || ca | cas 1 iD J ep etm Ih ip Vad Jams an a 

| | i } 


QIN NAIN © OUNAINN NAN HAIAICI 


ia 0 fae J Sn Y a a RS 2] | HY we a Ys Y | 

OOSGSOOOOOC OC eOGeia © 

ODDO COOHCONndDOO0Oe 

PIB LSMHKONS CUINAMaADHAMONEH O 

PR ARM MN DOAIA DO NAINA SLAIN SFU 

eeee#eetfeeteeee¢e¢#8 0@0 8@ @ @ 
1 


NANNANNNA AA 
O@eOeG lee 


Soeaeaeecse 
PD |G | SD I awk a 


NNAANAAIANANN 
COQOOoCcCOCo 


Seago coac 
ODVOOOCOO 
COVOONNODOO BDOCODONNODO 
WAM MAIR OIA AOAIMMINO©OO 
TAWA TOSANAAHONAMNMY FS OM 
@®e@e#ee8 ®©6UcMhC—SrhC<CNMMHMC—<C HC OhUC HDhLC<CMHhC<C MhC—<C<C MHDhC<C MCHC CUD 
OOVOODOVOTIVWOODOONIONONOO 
| 14 | 


MNIAQINIAIONIAIAIOI NANI NAIAIN etd 
OOO C OCOD OOOO OOOoO 


| ia Y Os of Ges Y Gay NY ee yy a Vf A Fam || i Va 
SOOO eaIOOOOWe oO 
Oe Sind Dera OO GOO 
FOWA SAHIN AM AOMOOANO 
m=O AUD OM AM FCS CHANIA O 
@®eeee#8t @e*#@ 8® @e 80 ¢ 8 @ ® 
6 eean | OB 1 |e ass i Ob 7 CF cae 9 a Ge ee a I EM i Fa BY ap 
| 14 1 } { Id 


NNNAIOQINAIAI ON Ne NAONAAN 
OO COO QOVUOCOCVOOOOCO®o 


Sr OOt a0 2 


eles a Ym cs f(s), a Meaty ae] a sae Ya 7 Oa Y ca fe (fn | [= 
QOG@OCUCTOOO0D OB Oeee 
OTVOOAMNODOOOOCCO90O000O 
NS AA ODIAM-LIAN TAG UORONAOD 
DDALPANUMAINALPASPALIM DOAAA 
e*eeeeeeeese#s#t © @ @6.O6hmUMhUCUc MhUC HhUlUCOhlUO 
OCOOQOOOCOOCOOe@Cee © 

| | | Lear 1 


NNOANNNANAINCIN AN AAINIAIA 
OO OVDOVDOOOOVDCODO0O000O 


SOooonoQooooomaeaojaj aoa e 
OO VAOOOODOAODIOMAOOV0V0O 
OCF OOOCOO2e ea ima © 
~-NOODTINAIAT SONOTCHD 
OMAR AN FiO OM OANA SIN OS 
ooereoe ee eer eee es ee eo @ 
! | | | i 


83 


LLUWS 


LU 


MATRIX EGS 


RCWS 924 COLUMNS 
f1CC 


© 


9 


MATRIX NUMBER 


OAM MMA SFO MH OAM TIAN OS OOH 
COOCOC8e CO Cea2OOorCc eS 


OoOeCeeeoacceeoace eae 
SCLOOMODZOCNOCACOCSCTOCMH!H 
aA OMMNORMONOTFOODRSS 
A SOMA FOO OOM AOQOAINAN 
TAT OSM PFHOn AKCN SAW 
e*eeeeee%eee?ee%e#eee¢e 80 ¢¢ @ 
oo rr ee 
| | | 


HAACA CAME  AINICACO COM OAM SOS 
OOCOWO0O CCOoGCeooqgeda 
CQooaCek. Ceoeceaeeeceagcee. 
OOMLSO CVWGCUODdVDODAMNIO§ 
OANIOSF DOH MADNS SFOKOMOM 
ANAM OA MNOUOMNACOUOOT DO 


WO NIC CAC OOM ath HS RSLANE 


oot eeoeoee eee oe eee 8 8 
OOCOCOCOG OOS ees 
1 | | 1 | 


TAA CACAS SFMOMAIM OM Meacaa st 
PFOOCECC'S COG eee ae 
OODOS24eocOOOeSsea sae oace 
KOMNMNODNODOCOCCI?: MOOOR 
DOOMKHPOMNOONNO Cmmninn 
HAHOSAONOAM GE MAMC NAO FOO 


eoscoeoveveveeeeereeee ee @ 
OOOOVUSGCeSOC ee eooeee 
iio l 


AAO COS FHM AG EOMNNNANS SS 
AOCC QO OOOO e ae 


AO MNOAS 2a Sera eaic) 
COT OUNVOOOS COG eOec Oo 
AaNeItPNOCOOVNOCOCNTOWM 
Mm OAM FOSNAOMLODNMAYW 
OO aN TO et FANG SO of SR 
e® eee 8@8efeeeeeseee#eseeseese e 
GOogegeac tO SS game 
| 1 | 


LAVA AY CONS Fs: 
OCGeoocoeg © 


OY EN T GVENI TON] 
OOO ae. G? 


SOOSeG2 Se] SOs eaea 
OCOD ac OOCGC Ce vee 
NAMOFHROWO NOCDONOW 
TUITIARAINOMN NOMOONAO 


MA SFORNOOMON CNT NeMON 

ee @ © e@ ®@# @e8t# @ @ @ @e@hUhhmUCUhUhHhmhUC<CNMHhC—~CM Hh—~C OhUhlU!F 

CCOCGeGi OO OSS ee eee Oe, 
| bal | 


CNCO SE FNC ONO OY ONE VON OA 
OOOO VOMSOCOCCOOO 


ee. 


Soa. waSoO ao eao Oo 
TO OKO VOOO COOCONC @ 
TOS CIOOUNESN STAC TOOMNN 
TOM KHON WHAONMO FA 
eeceeeeetee ee ee © © @© 6 
COCOC Ome SOOO OO Be aes 
| laa | ! i 1 | 


VERVE ON APE A COO EO AICO Oe OO 
YWUUCOOCOCCOOC OC OCO00 


SCAaCCeoGeaSCS ez eaeee aS 

2OSTONNDC COO UOe2eCeoea2o® 

DAT OT OWNOBDNDDNO SHH 

NRO ADM SMM BOK ST SOI trooe 

DM Dt Da MSIN SENSE OOD 

eoeetee»ee##e#eeeeeeeeeeeeeeee @ 
i | I | 


CNA ANF OMA OACAMCACA AIC Fw 
| aie T ibe 3 ap Yam Ii eX iby Ye i ee | SV al ap > ew ea (eo se ee 


ODOoodOCooCcea seo eeeasSs 
OOOO OVO DVOOOOCVer™ 
PN OON M FT OWMNW DO DMNOOCMNON 
h DUNO DMHWMOAT OOAOOMWND ON 
TOON SD NIA SO at 0) et et et 
eeeeoVeseveeeeeetee @ @ 
OODOOOO OC OO9COCO DO OCe 
| ] | | ‘ | 


0 eJLO000 
C.G1L000 


TN BLANCA = 
ERROR IA BLANCA 


ERROR 


co 


4 


MATRIX U FOLLOWS 


6 ROWS,+24 COLUMNS 


4s 


MATRIX NUMBER 


NNN N NANI NIOUNUINN AI RIA OQU AL OY 
COOCCOCOCOCOOCOCoOeCCC ese 


SCoecnonoooaoooooaneaak. 
MANOMMOAT OMIM POT C hm Sete 
NAM AM FAW O ONO OMNANININ 
CAM SAM SANDE OEE MNO 
IPL FINA FOOSE OPM ON SBP 
ees © © eee ¢€ 6 4 8 6 Oe eee 
(oe ae 

| | 1 | 


ONNNNNONANANANANANO OU 
OOCOOCCOOCOOOCO00O 


SCcogcneaoGooaoCcooceeca 
OH OMHAOPODMATMNAOONO sO 
AONLKOCTTOMOORHINM af 
OR ONOMOMAODAANELE AO 
MAIONITE MON SN PCM HN CONUS 
eee#e#eeee¢eee9ee#ee¢ee¢e¢e¢6¢ 6 @ 
eT Oe Oe OS aaa 
| il 


NNAINNAN NAINNNC NAN AIA 
OCOOOC OCOCOO°OCCSCeae 


SoOooOOocCocoqgcKkoecea 
Set It DO ONE AR OUND ONO DMALAN 
CO RIOR STOOL ON 4.010 + ONO ea 
mt FO HMO SACHA ROMMAMINAC 
ee et Oo ee 
e ®*.6h6€°é eo ¢e¢ @ @®eeeee® ee e e 
ODCOCOODDCOCOCCOCOCS 
| { ( ia | 


NNAUNNONNNNNN SINAN 


BOOADAOCAOOOOOOeaa es 
SLINF tO OD NIN Sam OAM ARAM 
=SrHtOStTNAODAMNSTHOTOOM 
CAT AM LAL PS LA TNO HOM Se & 
ACAD OPFOR ANNAN ANNA YS LC 
ee @eee 6eeeeeeee#ee#eeee#seese ese @ 
SO > a 
\ | | | 


NNANNNNN NAAN HO 
ODOVOOCOCOOODOCOOOCCgse 
| 
COOCTC SAOAQ0OOCOCeSa es 
NEEM AMONNE TORS COMA 
ONOMH MF OGOSMOMANN CO 
OP RDOANARSOOPOMAN aM se 
TPOP AD AN AINAOM SOOM 
e® @ eeeeeeeeee#%eseeeeeee @ 
CCCOODOCOOCOCOOCOGGae 

| 1 | | 


NONNONANN ANNAN OOS eee 
SOOGCVUGSGOCCOOGCOOOStCCO 


SBOODOOO0O S00 Cceoeeas 2 

DORANAOMT ONDE OR SIAM 
CAM OIN NO A OO SUN OL Daa 
COM BOG TD OOM 01.0 Ge 
NIA Neat 0 OOO OG ae OSL NIP OO 
e* e@e@ @ @® @e@eeeeeeeseee#tee @ 8@® @ @ 
ClCOCOCOODC VCC C2 ae 
| ' | | | | , | 


OOVOCOCOCCOCOC COOGe 


NADOODOQC203D020C6G8ge 
CrYPPOVOLINES thar aorsr 
NOH OANDNODNATARN™ oO 
DO AFODDEMANAHORINMNACE 
DROS NAN FON Dean 
@eeeeee*eeeeeeeeeeeses# eee ee 
QOCODDOOO0COCOOCOGees 

| | | oe | 


ANOANNNNNNNNNRANNAIAN 
ODVOVOAOD OD OOOO OD0OO000 


Ooo0ocdooqoo0008steae 
MIS OP TOWWMNAINOR CAM OO wr 
LN OAL mt SION CO AIR BL ST OSE LY mt 
OD ADUAMMOARASONDOAMO®D 
MAINILN TD PUNCVUN DO eat O AN MAP UNO SF 
eoeeeeoee#e#eeeeeeeee ee ee @ 
OOOO VCOOOOOOCOOOCOOCe@Gese 
| | i | ; | 


—_ 


LU-SQ@UARE ERROR 


2: MIMIMEZES 


CUTPUT NC. 
GESTREC MATRIX IS 


SU CCLUMNS 


6 RCWS;, 


EDs 


MATRIX NUMBER 


ONO) ad ot et C\) 
genes 
| 
CeEeceEede 
OM TOO 
CEA DH ACN CO 
Cuca for 
Ne vatvate mat 
& e.6.h6©° e.6h6©°8¢ e 
CCeBe eG 
| ' | 


CORFOMeH 
CC@aaG G 
| | 
CoB eae 
Troe com D 
=F ORG 
CINTA OM 
wot oT OUM of 
e e e 6 6 e 
COO Gea e, 
tt , 


CAaACNOr 
Oc wee 
| 

COB as 
Ci Ola = 
MP AM ON 
AS NSO 
CO) mt a eet NL 
ee®*e ee 
CeCe ae] 
| } | 


oad ed md Cd ed 
CIO OGa a 
1 | 1 | 
ais ss Yb aan la) 
NT Stim a 
Nya 
LN Oy FN TE 
nO le CORN ot Dera) 
oee eee 
OMaGte = 
| { 


ONaAOOND 
Coe ae'c> 
t 
(JRA ES 
SP DM wT fh 
Qe eoens 
COO Ot N 
tO) et et 
eee *® 
LB GB | To | a ew 
| | | 


aADHOOD 
VUMOCOO 


| 
a es js Vi fa 
mM -rasienc O 
WN Oem A 
NONADO 
Cet AOAAS 
ee ee 6 @ 
E> ee es Ge ee, 

) | 


WS. 


See 20530 


y 
\ 


EQUALS 
2 EQUALS 
3 EQUALS 


NUMBER l 


CRITERION 


2 ow keaibZens 0 3 


CRITERICN NUMBER 


C6o7654D 04 


“2 
— © 


CRITERICN NUMBER 


ELAPSED TIME 


] «G5 
MINIMIZES TRIANGLE ERROR 


DESTRED MATRIX IS 


85 


CUTPU NO.3A, 


6 RGWS- 6 COLUMNS 


L224, 


MATRIX NUMBER 


=—aOODNN 
dean 
Saaeaae 
MA $OhMAIM 
P= OO) > &4 CO 
GOAN O MOO 
AI set ef CQ! 
ee eee? @ 
OGaGee ee) 
| tt 


CORON 
SO a BGS X so | ib} ep 


| { 
OOo aes 
MO omOmon 
MOD ORD 
AIAN tA OD 
a aONS 
e® © @h6mehlUD 6 
COCVWeC eG 

1 | | 


OOCONOO 
OeoOoooo 


Ge fsa | ae | a | am Fh aay 
Veo he SEN ake 1 tag 
Urs Or OS 
af D ea M 
COA = CO) = 

ee @ @ @ @ 
(> GD | GB 1 Se SY 

' to 


aeOaBMOAO 
OOO ea a] 
l ! 
LC jf eas (a oa foam fe 
<T oO UN CO Ph CO 
CNNO aM 
ees e@ 6 ® 
COO OOO GC 
| t 


ONODOOO 
OO0OCOO 


ACAD 
INeSOOS 
MON SOM 
ODN 
i ed aN) ed ad 
ee @ @8¢@ @ @ 
Soe Oc 
| ' | 


sqaOmnOOr- 
ODOVOO 


| t 
OO0O9008O 
NFTNom 
SOON 
MO NM NO 
DANA 
eeeee?@e 
OOO9d00 

tt 


°©1454270 O05 
Z@192255D 


a 
=) 


EQUALS 
2 EQUALS 
3 EQUALS 


1 


CRITERIGN NUMBER 


G3 


Oo} O04 


NUMBER 


CRITERICN 
CRITERICN NUMBER 
ELAPSEC TIME 


™ 


2 ek 


1.61 


MATRIX V FULLOWS 


G 2UWSs24 COBUMIS 


MATRIX NUMBER 


OA COM AHI OMA COMMAS OM 
CO@@COC OOGCOCOOeCee eo 


SeaQggceaceoounoaeaoge ase G 
NOCOMATNRANOCMAOMOMN 
NOL HAF DOMACORANME AN ih 
FMR MDS DA TN DOAN N et 
FINA SFHAPAMS OHRON SAMOS 
®*e0e08t &©@h8H 8 & © OO HhmUcHOMmhUC<CNMHhC MHOMhC—~C HHMhC—<C M HHhC(CUhHhUCUchmhUllUhF 
ie Gn iene te 
I { { 


ANON OMG CAIN AMM MMM SOOM 
OOMOCOCOOCO OS OeSeeee 


Hocoackcoceceaqeaagaac 
COLA mT OD OS ORHAN ONS SOM 
MAR OQNNANCOOSBOOF HEE OWM 
SNA FOM-OnR OE OaCDOONSTE 
CSANANAOAD Oat OM Sat AMm O&O 
e*@e%e*ee?8e#8ee08806hHGhCé<‘<“°€Lca¢ « &* ©~« eh ehUh!hlhlUhO 
SS eae gaa 
| i | 1 | 


ORANG OSPHMO OI OAF OMA SO + 
COOCOC CO CeeeeeGeeoe © 


Aaaoogoqgagaqceeceaaaaeaaae& 
= NOM AMA LTAMCOOMAD 
COO SRHAIN OM AOD HHAIOm ONO 
OO CM SCOR Em OILT CLA NS CD 
Net Set et IN OAILN FOF SAO HM SI 
eetee?ete0ef @e@¢©6.?*®6Uc!8OMmUCMOrmhC MHhC—<CNYHhC<“ HMC HMC UCU 
QOOOC COG a ee a1 OOO eS 

| | bt | 


AAAS FSO OAM STS MOANA SS 
COCO OOCGw es OOOO oOo © 


QOOFPCa2 oS Can ele eae Ss 
OMmNOEENDNOCIADMAOOS™m N 
Ore OOUMFONNATEHAOOONSTH 
UN TOME MANAIE OOM Ost Os HO 
OD Dt OO tT IAIN FO NF SI 
eeoetee@ef#@®@6mUh/Omh6UhOO OF 6 8emhUhMHhUhCUhOmhUCUhHhmhUC< HClUh}}HhhlhlhO 
OS a eons tamer 
! ) | 1 | 


PA COONS OA TOAIMALC MOM OA SOO ST 
OODODOQDOO COO COO OGOo) 


QA CvSae OO ae eee] 
PHAISP TA CEN AI OOK ANO- 
TM HM OMA HON OAOINO © CO 
OMONO BOM PFOON ALAM MNO ® 
OS PANO NE DOAUENONAN 
eeekeeeeeee8« e®6UmcahcmhMhCUmCUrhhktlmchlUlUcrOhlUCUchlCc HhLhUlUrh/!hlUwO 
COOC COGGeaOsDe@QQaoce]) 
! I ! 


PY OM MTT MAIMN MINE EEN EO 
GS LE I ay 9 ET Es TS ee fk BY ie 1 OD Tt ape Yee I ee 7 a | 


OoSSeS a C2SOo Og CC Oeee cu 
TV DOM TH DK AMIDA N OD ONO SN 
DD Cet ON Se) SO SOL OS CO 50 Core 0 
LEN SO Sat Na = NON CCN 
ht tt ST RIO mt ot 
eeeteeeeteeeee#e#¢e#ee¢ @¢e e 
O° COO ea Sea aN OO mes 
' I ! | ' | 


eaten ton tan loataa tea tanlabtaalaaletlaalaal laakaaked: 
QOOVOMDOVOOCOCOOUYDOO0O 


Ono eas{| Se aeaSaas 
AVN TIO SAAD TOP PR IO tO OS 
et (INP OONO MD Aa OD DNO 
ODARDNOPANDA HA ANAS OME Mm YF 
CSDMLPM PO oes rvuSPOoMmwnc 
eoeeee8. efeee @8e@©6hmUmMGmUCUCOChrhUC HhUCc HhUCUc OhUlUhhUO 
OC0°0! == ay eet 

I { | 1 1 | ! 


NOV SAO SP PFOMADMON AAO ANOS FT 
OOQOOOCOOC DO OVC OOCCCeeads 


SD AICO O Oooo ose SS 
WOTROML WIM OH UWUANIAT—E UN 
DANA N SA CHD HCMOY SALINAS 
AT ANDOOM ANAM Pet NWDOA 
DAS SD wt Ht NI AIO eA tN tO) St et 4 
@eeeeee#eeee?eseeee?se#* ee © @ @® @ @ 
GSOoDdoodQO00OO COO eGeeo =] 
| I | 1} ! l 


86 


Lee eros 


ELAPSEO TIME 


QUTPUT NC. 
DESTREC MATRIX 


LyMINIMIZES V-SQUARE ERRCR 


TS 


6 COLUMNS 


? 


MATRIX NUMEER 


QOOO AN 


| 

Ceada eG 
SON Oe 
On Oth 
BAU OO OS 
ee ¢ @ @ @ 
3] at ER ee KK | 
| | 


Ooo 
O@GeQgac 
| ' 
ae eee) 
ce OMr~ro 
Wan ol mal Of dl De 
IN 2 PFOWCSY 
—a HAM eS 
e® e0«© © @ @ 6@ 
Sew o©® 
| 


OOCNODO 
COO Oeo® 


ei a ES af Bf 
rrm9norr 
Owosremoas 
WECO DN at OMT) 
OT OI aN 

e @ e e eo e 
Oooo ®> 

| 1} 


saOHxIO44O 
OeOooU 
] I 

is Fs f Cs J as Jf Ba) =) 
OE Corn) 
TT SF Sirs 
CSOs = 
teat PN OY 
ee @t¢e 8 @ 
Qe Co Coo} 
| ! 


ONOYOO 
COO000O 


SSsaQgogat 
pA PP 
On wsosrnr 
MOR DOAON 
at md mt (OO) mt 
@ eo 8 ® @ 
Le Fey 1 Yi 5 ae 
f { ! 


geet San SI GD If | 
OOOO0OO 


! 
as 1B a $f as) ef 
MMO TON 
0D DAOWO 
NIG NUN 
Oar kar | 
e® «© © @ @ @ 
ODOD000O 
t | ! 


4.9981920 05 
222846200503 


1 EQUALS 


CRITERION NUMBER 


EQUALS 


2 


CRITERICN NUMBER 
CRITERION NUMBER 


ELAPSED TIME 


1459580 04 


2 
—_ @ 


EQUALS 


5 


1.47 


MINIMIZES TRIANGLE ERROR 


OUTPUT NO.3E; 
DESTIREC MATRIX IS 


214, 6 ROWS, 6 COLUMNS 


NUMBER 


MATRIX 


ADOOANN 
OO we 
| { 
SOQaeaQaO 
CA S$ COM AIC 
P= CACO FF eID 
DO ALT AD 00 0G 
OS eed ot a PS 
e® e0« @ ese © @ 
\ tf 


OOS a 
OOO COG 


| | 
a | Ep asf | 
OOM OCA! 
CODD N 
AICO ALA CY 
aa ony -- 
ee *%* eee 
ee as 

{ | 


COONS 
QOGBo 


Wise CO<f a 
LOU SOE 
matt OD et ee 
COV ONS med mnt ON med 
oe 8 @ @ @ 
BY, 0 EG CS OS 

N50) 


AO et DH CO 
GF 1 Gp f GBs See aE 
y \ 

i etn fa | Se 
WT OLN Oh 
ONIN oon 
Nam DOL 
CO AS HO -+4 
® e °® e 8 ®@ 
Ot Cate 
{ | 


LOT ga heh Ep | CO ae 
GE GEN Sele Yh > i Gk 


OSC 
TNO OO. 
NONSFOM 
DGOHmaAtPOal 
rad md med CAS md ed 
s e t ® ® © 
OOOGe 
/ { { 


ADAOCH 
Ogooudo® 


| { 
Nooo Da 
Nrrmorm 
m=) OW Dh 
MONAN O 
CAO MAA! 
e®@¢@8?@C=«UmM—=~<C CD 
OOOO Ce 

1 


O5 


145427) 


IS2 
2.0127C5D 04 


c 
~ @ 


EQUALS 


NUMBER l 


CRITERION 


ee 


Ze 2590 


EQUALS 
3 EQUALS 


2 


CRITERICN NUMBER 


CRITERION NUMBER 
ELAPSED TIME 


Ga.0 
MINIMIZES TRIANGLE ERROR 


CESTREC MATRIX IS 


87 


OUTPUT NC.3C, 


Zils 6 ROWS, 6 COLUMNS 


NUMBER 


MATRIX 


MoMODOONN 
i iaeaaeaiua tual 
ec fs 9 os ff a 
COS COM NO 
P= COO DP ed CO 
SO ALEVE OY 
NI ad med med P= ONY 
ee ee @e0 & @ 
Oe oe) 
| 1 | 


8604.0 Go GP IT oe Ne 


7 | 
sip | A in 4 a 
TO MOA 
DLlOKNKNRH 
QE COLA at LA OO 
as Oam hy 
eet eee 
Se ais 

14 


Ge Ne © 
ST ie Ma YG ep ey 


(22 A aS i BY 
IN 0 CO wf OH 
Ls SP SOM OS 
at Ons aeS 
NAS aS eed A] OS 
oe#s#see® 
COCR © 

i Vf 


at Oat OrtO 
OOVUOO00O 
| | 

Soooaact 
POUND © 
CIONNOAOM 
ON at TOLA LA 
OAS HO 
«ef @ @ © 
OOCOR©O 
| ! 


WN OGamO©O 
OOO Om © 


Jy a i YY 
ITNDOOST 
MACOAIST OCS 
OD ASTON 
aaa A at 
e® 2° @ es 8 9 
a Ley A 1 ED ND | a 
| | { 


AOMOO- 
OdO000 

{ | 
Ooo G@as 
AIST STUY Oe 
eA OO Of 
OO ONAN O 
Dat OAOAAS 
e° @@8¢@ @ @ @ 
Oooo00o 

1 f 


521454270 05 
Celio 25 208 
42.C127C05D 04 


1 EQUALS 
2 EQUALS 


CRITERION NUMBER 


CRITERICA 


NUMBER 


EQUALS 


2 


CRITERIQN NUMBER 


MINIMIZES “TRANCE ERROe 
IS 


OUTPLT NO.3C, 
CESIRGEL MatR ix 


é COLUMNS 


6 ROWS, 


253% 


MATRIX NUMBER 


=OCONN 
COOCa® 
t | 
CaoCcCaa 
ATOMNGD 
Ph OO) cf wat CO 
WO AULA O09 OS 
QC tat et A) 
eo ee @ @ @ 
COCO SoOe 
| | 


OO aa) 
OOIaaco 

| f 
BS ib Op | CS A eB 
C2 os ON 
TO OO OO 
AICO LD tL OD) 
eat ad SNS 
eo ee @h6l 8 Ul 
eat sa 


QOOON@ © 
OO@CwOO 


2 GaicnGs CS 
UNO Oxf OM 
LAS SOM OS 
mat Oma) 
CO NY et ed NY ot 
eo © @#e ee @ 
ODO0CM%O 
| ge 


Ont OmO 
OOOO00 
| | 

Oaad000 
Wo WV co & 0 
ONNO OM 
Na COLIN 
OO) med ST met OD ome 
oo © @ ee ® 
OODOOO 
' y 


CONOOODY 
OOOGOD 


(a Gary Gia) A | am) 
ITNDOOT 
MNOAIST OO) 
O Dra OOA 
ed ed eed CN] md md 
ee © © 8® @ @ 
OCO@oo@ 
f | | 


aADAOO— 
8 | Ge ES) Sb | >. 


t 
Soe = 
NTTNOM 
m= OWN Om 
PYVO Net SO 
DAA MAA 
eo @ @ @ @ 
OOTO00 

bo 


59-145427D 05 
Cel Go 5> Das03 


1 EQUALS 
2 EQUALS 
3 EQUALS 


NUMBER 


CRITERIGN NUMBER 


CRITERICN 


88 


O1lzg7C5D 04 


2 
~ & 


CRITERION NUMBER 


ELAPSEG TIME 


109 


= 18.Cé 


TOTAL ELAPSEU TIME 


C-2 Sample Problem Usins Progran CECI 2 with 0% and 10% Maximum Random 


Errors. 


The nomenclature used in the output of the program is as follows: 


MODE=1 


BETA 0. 


Matrix 
Matrix 
Matrix 
Matrix 
Matrix 


Output 


8 


ig 


BE 


No. 


Data will be generated in the program 
(The sequence used is Ne Ss Lee. Oe) 
Original matrix C as read from cards. 
Matrix U as read from cards. 

Matrix V froa V = CU. 

Matrix B& with random errors introduced. 
Matrix EB with random errors introduced. 


The desired compliance matrix which satisfies 
criterion four. 


The printed output for the problem with zero error appears on 


pages 90 to 


inclusive; that for the problem with 10% maximum 


error appears on pages 4 to 97 inclusive. 


89 


MATRIX P FOLLOWS 


6 COLUMNS 


6 RCKS, 


ls 


MATRIX NUMBER 


OOOD Ol 


OOO000 © 


OOO00 © 


200D 02 


OOOHOO 


OODOO0O 


O0O0D Cl 


OOWMOOO 


OOVOOOO 


MATRIX @E FCLLOKS 


~~ 


6 ROWS,24 COLUMNS 
-C.19CCD C2 


29 


MATRIX NUMBER 


NAN Ne NNN NN SIONS 
elelejelelolelelelelolelolelelelelea 


AQoaQoacoaeosaoecencooeca> 
OCOOOVOCCOCTOOOCOOCC°CO 
NPTOHOHBOOMOOAOODBDOOW 
ONAN IAE NOD AUN SF OUD 
VO SPIN COR OM OO CON ett od 
eooeeeeveeveere eevee ee @ 
| ' { 


NEN NANNNAINA 
OCOODOO900000O 


OOCOCGOaG0g808Re8 
elolelelelelelejlelelea 
COoOrcow DWOOODdOCcoON 
ODOOMS NNAANO MOO EC 
MANO NSO OOO COMO ONUM 
eooeoeoev eevee ee © @ © @ @ 


290 OS 22S mae | ~ atm 
rt { 


ONNNAIN 
OoOOceo 


COOQQ000 
OOO000 


NAINA 
OOOO 0O 


ooaceaS 
OO000O 


NNNNNNN NAO 
COOcTOOCOCOCO 


Ooaoaooooaoaos 
OOOVOOO9OOCO0O 
MOMNNnOAMAOOODCOO DCOOD 
TNOO-MMMODOHOWM WOME 
Nam NY SOM MO TIN ME ODM MO ae 
eoeoeeoeoeeeeve ee eee © @ @ 

ro { 


NANNNNONNNNNN SA ae AAA 
OO COV DOOOO OOOO 0000 


Sooo ooOoeceaoaaocgeeqaoedce 
elejlelelelelelelelelelelelelelelel® 
OO DODO VOHDONMOO0O00% 
BATH KHON SOW MetOet Muh O 
AR AAPM LA OANA MAN DN LAU 
oeoeeoeceoeee eee eee & © @ 
! ' 


NNAINN Nata 
OOCOWONO00O 


OOoeago0e600 
OOOO0O000O 


NN AEINAINAIAA 
OO OVOO000O 


SoQaeaQgoks 
OVDQDO00000O 
OVOOOOIMOOODO WVDOOMMWOO 
Wee AM AOA OAM OO 
AINA HN HHO NIN FOLN OOH SF OO 
eeee?%eee#e@e@6mthmUhhmUCc OhhlUc OHOhUC DC‘ UC MULUCUCOhUr;!hlUO 
ODOODODONO OO DOD VOD0O00000 
in | | 


NNNAINAINAAN AION AN A iN 
OOOO VO VOAVOMOOO000 


OOOOO0O00 00000000 
ODOWDOVOVOVOVOOO 0000 
ODONODOODOVOVCOO0O0 00 
=OWIA ST SANA HDMOUIN Os 
OPM NAO Mat Met TOO MIMI AIO 
oeeoeeeeveee eee © & @ @ 
OOVOPDOVOVVOVONOCO00O 
1 ! t | ti 


NAN AIAN AIAN NN eA ON AIAN 
ODODOOO OO9O00 OOOVO000 


OQOOOQ00Q00000Mm000000 

COODOO9O000 OOQO9O0O000O 

OT DOKRRNODODOVOODOOO000O 

NOAMNOWIM IAT HAODRDON ato 

KNOAPAMANNMOA SOS NIM OMATA 

eoeoeeoeeeveeveee ee © © © © @ 
' rt 


NNOENNNAINAINN Net NAINA 
OOV DOV VOOVOVOO0O 9000 


OQOoOOOo0QO00000000000 
COOVOWVWDOV OO OVVQDWOO0O 
OOVOFF OOVODVOCOONMOO 
~-NOODTNHATTOMOTOMD 
AANA AW TNO ADOAAIN TINO 
ooeevo eevee ee ee © @ @ @ 
ODOO9D9O OOOO DQOO9O0000O 
J ' | | | | 


90 


MATRIX EB FCULCWS 


6&6 RCwSs24 COLUMNS 
-0.171CD 03 


3 


MATRIX NUMBER 


AMAA GHNI SEMA MO MOMMAIM FAM 


CeaSaSaceoooooceoececoesS 
WMWOODDODOONMNOCOWMOOTOON 
ate AMOMMONMONOC PTOOR STO 
WODMTMMDSMOOMOOP HDOOANAAN 
FUVLHR Pass HqMnreQuoneqaniw 
@eeoee%rk#ee?ee#e?@e%eeeoe#?e%8e@¢e @ 
ee ee eee 
{ | 


NNOAOMA OMS OS 
OOOO COOOCO0O 


COOCCO0CAB0O0O 

OCTOOVCCOMOrs 

ONOFT OCH COBNATOLOMM 

MOUMMON NODMDOOOST DO 

DANES OD 9S Oe met oe et et LO ed 

eevee ow eo eee eoe eee 0 © @ 
{ { 


MAIM OS 
OOOO0O 


COoocea 
OOOOer 


AOA OM 
OGoOoO0o 


COoeoceasa 
OomsToo 


KHOAO CO FLFOMAUMAMG MH 
Oooccooccoceoq 


CACOCOQOCB0QAdDO 
MmOWMODBNAGDOOOO 
DMOaHROMMOONNO DOONnNnN 
=HOACNOMOMAHAO NDOLPOO 
QUO et St HOA O SF OO HOO Sf ot 
@#eeeeene%seeee8eee@¢@ee8ee8 8@ @ 
eo 00 OOO aaa 

{ | | ! 1 


MOO SPH OMAAFMOOAIAICIS ot 
elelelolelelelalelelelelelelelelele) 


Sooocooeoeceooeoeezoeceas 
OOTONOCOHMOQOOVORALO 
AN FNOOCONODOOLOOM 
AMO TOaHNAOMPNOCNNAY 
COD wed at eed OS OY et SPO OSS OO SP met tot 
oe eeeeeete?eee8% 80 &©@ @ @ @ 


OONPOMAIOANL AMMOMATM + 
OOOOOVOOOGO WOOCO0O0O0O 
OxqooogCooCese aceoeaeacke 
OODOCGOOCCOO DWDOOCOCOWON 
NA-OSTRAOMNO NOVODNON 
TUDMINMAAAOCIA WRWNOGNAO 


APO HAIR AAIM ONON TVA 

eoetreevevereoer eevee ee ee @ 

OOOO OMDOOVDOVDOGVOOCOC00 
' 11 1 


NAF PFOANUAOAMN MAN OMMOMOS 
GDOOVOIODOMODOOVO0O00°O 


Hoe eqcoeooaeesaaagdsasS 
TOO ChOCWOWDDOWVOMAOANGE © 
TOHMOOMONTNSFOOMNN 
TOM HNN DHOWOT AO 
DO 00 ed ELD at mF COLD at SO ett 
eoees#ee#@?8#e%@®@@?e?ee2e8 @ @ 
8 Oe | Oo eee 
i \ | | { 


MNOAOAMAMNOAOA AOA NO Oa ACS 
OODODOVODVDOOVOGCVWOCO0O0O0 


COCOOGOO0O0O00000000 
ODOTSTONNOOQDOCOVOOC0O0O 
DATOTOWNODNODNOSLH OS 
NA OAOTOM HOO ST AO OO 
OM TOT QS STIENTMOMODYE 
eeeeet@eoeeee8e © @ &@© 8@ 8 @ @ @ 
2 OO ES eS See 


AMAOD LFOOAMOAMMOAMAIMS SL 
OOO OO VDOOCOVDVVMODOO000O 


an] an] an] oY onl an] oO] en] en en] en] a nl el ae) 
OOOO DOOODVOCOOOVOAM™ 
MOONEH TOMNANDADONOOMANON 
PR UN OR Om LASS OO DO OWN OM 
ONY Fd et DN) ed Be eS eet OF) eed raed eed 
eooereevoeeoevoeeeee © 8 © @ @ 
| | I 1 | 


91 


0.0 


ERROR IN BLANCA 
ERROR IN BLANCA 


0.0 


— 
= 


MATRIX U FOLLOWS 


6 ROhWS,24 COLUMNS 


41 


MATRIX NUMBER 


SNANNNEBAINNINAIANIN Ae AIAININ 
C$TOGTOOCOOOOCOCOOCOO0OO 


Sa cooOoeoocooocooecoao 
—laololololelslolololelelsiololel> 
NFOHOHOOMOOCAOODOON 
ONMNNANENODAMINMIALOND 
INS TINNDOMOR ME ONAN 
@e*eoeserteeeeeeee?eee#eeese§e?e?ee@ 
FSOGS0099S00990005 

| { | 1 | 


NeAAINAIOIQUNIAIN 
COOC9O900CeO°0o 


anf cam fam Je cam on cam foam ei fa Ya 

QOOOOGCCO0O00 

OCOOrDOW DOCDODWOOROON 

ODCOM LT NNO OHO 

ENON SHO SRBMMHODAENOTAUM 

eecereoeeeevee eee © @© © @ © 
{ | 


NAINIAIQINININ A AIAIAT 
OOOCoOoo000CoOO 


SoZ Soooecoeocke& 
COODODOCVCOOCVO00O 
MOOWNKOHOOCOCDOO DOOVOM® 
TAMOMMWNOCHOMW NWONWNWP 
Nat At OPO OPIN MNO SOM HHO 
@®eeeee??etgee¢e¢ese@e@6mvh)eee8e @ @ 
a 
1 


NAINANINONAIAICI AINA AIEAILCIA 
OODODOOOOCOCOOVCOC9C0000 


Oooo condooononeooeaca& 
COODDOODVDOVDDVVOVOVOVDOOC0OO 
DOOODCTWODVOHVDOONnNOCVOONOOM 
AAS AOWNSTTONGOHOHaMONlO 
DM AHO MINOURA DO AAD SUVA FU 
®#eeseeeer?e@e 208806 8@ @ @ 8 @ 
— oT Tee 
it if [ 


ONININIAIGI 
OCcoOo0o0o 


COOCcaOo 
OOO00O 


Q2 et AIAY 
Oo°Ccodc 


eo S 
Oooodo 


AQINAAIOIONIOAIOIOL NOAIAQINQUAIAN SA 
OOOCOOOOO DOOOVO0O 
COCOCOCAOdaB AGOeOoododoo 
COQWOOCOOCOO WDOOOOOCO0O 
OOCOOMOOO DOOONMNGO 
OHM AM AARON = DAAMUOO 


OUD SFO HA HLA ROCA COLA SP COR 

ee eeeee?e#e?ee0e08808f08086@t6hmUmOMmmUMUOCU® 

OOVOOD OOO OOOOOCO0O00 
I i 1 { 


MIAIQIQIQIAQIAIOAIAIOAIOCQIOIO A tN 
ODVDODVO VA OVOOOOO0 


OOooCoo Coo ooeoees 
COOVDODVDOVDOVOVDODVOVOO 
ODDOONMNDVODOODVUVOODIOVOOY 
SOWA ST Saal ROWUINOS 
ad OPM AIA OM at Meat Ff OO SLAANIA O 
eet eeereeee eee e& © © @ @ 
VUDDODOVDVIOVSGCCVOCCOCCOCO 
{ 1 { | ! | 


NAINIAINIQINININIAIN SON ONAIN 
VO DOGDOCVO VOODOO OOCVO00O 


os. OCoOOOOOOOeaeaaeco 
OCOCOCOGDOCOSCOOOCOO0O 
OTOORBOOOCCOOVCCOCOO 
NDA OUE UN SF HOODDON AH 
AOU SF USM AUM FO FAL MAO MAIC 
eoeoeoeeowreoe eee eo ee eo @ 

l 1 itt | 


NN OBNAIAQININAIAIOCIAI RA OUAILAICY 
OOOO ODODVDOOVOOOCOOO0O00 


CORBDDOAOONONOOO000000 
VOVOODOCOOOVGOOOO0C00 
COOOFPOOOOCOCOOONMNOO 
MmNOOODPNAHATT CWO TOMD 
POAAILPN Oat OMAN SUNOS 
ceoeoeoereoer eo eo ere ee eee @ 
elelelelelelelelelele!s)a=lelelele}@ 
| | { {1 | | 


9COQD C2 


MATKIX V FOLLOWS 


€é RCS ,24 CCLUMNS 
-Cel1L/71LCC 


Dy 


MATRIX NUMBER 


HAAN MAL TOO MAM MOM NO SOO 
OOCCoooCcodo0cococeood 


QOoOQoOOnOAeAOebeeeaeaeaae 
MWOOCWMCOCCMNOC@MOCcCYTOOY. 
SHAN OMMOHONOTOCA SO 
WBDOOMAIDAST OOOH MOM BHOOMINA AN 
PINS OS mt THOARMOAON HNN 
@®estesesee 8 @ @—6@6hCMOC—~CN—é~<‘ HUC~CHtC<C< HCUMSHhUCU /!hhlUwO 
habe Ga ee 
\ 


AMOAMMOM ANMOAMMOAMMMST OS 
OCOoOOCGCO dOoOQcodoceQqnom 


Seaocca cooockceeeezeee 

COOMrFTOO ODCOCOCCOAMO-4 

ONOFSOSH GTODANSTOLTORMM 

TMMON NWNODODMOCOOFDO 

WAI CE OY OOP ett Em Bete 

eoeeoeveeve eee ee © © © @ © 
( 1 | 


AMOAOAST SFOAMAMMOM MNMOMOWS 
OOMOONTOOOCOOG C00CO 


EGoacooooooaoceaweea@ac 
MONOLCNOOOQGOO0O AX20C0- 
DOMKOCNNOOAING AOdMNonn 
AOAC MMHOMeIMO NOSTCO 
NAN See SaaNOTOOROOS A 
eoeeee#e¢e?e%e%%@##ee%e e# @ @ 
oe ee aad 

( \ ‘ 


MOOS FAMMNMOS OMAN SS 
ODCOOODMTOBDVOOCOOOCOO 


2QaCcCcOonoannoooeceaeega. 
OCytT ONCOOMNOCCeCeeeoaG 2 
“MASA DOODOONOCDMOTITO0M 
CM OM™ FORANROMNODAINANY 
OOS AAA HASH AIAISS Of Hee 
e@®eet%t%*ee8 @e?08@8@he@h6UmO8mUMmUCUOHmUCUMOmUCUMSUCUC 
ee 1 1 ee 
i | | 


OAS OCONOAL OCOMMO FOS 
COCCOOO00O0 OOeQCceoso 


Seccoaoaeea CcCozeezeaz] 
COMOCCOO0O0O DBDVMOCOMON 
NOMOT ROMO NCOCONON 
TINTMNMAWMOWN NOMOONAO 
OF OMAID AAUP O ADO NST NANA 
ee © eee @ ee @ @ © 6&8 hh HhUhHh(C<CMCUr!lhlhlUO 
ODDdODOCOQDWOOOO00OO000 
' | t 


ON OAS SS OV ENCE) C0 Os ae Co se 
OOOCUOOVUODVOCVOO0N 0) 


‘Oe 


OOoOCOoo OOo cooOe eS = 
TOs DWODOCOGOOONODD 
TOR QOONGNYNTOOMN 
tOM AMINE City =e 
eeoeeeeeeeeeee eo © @ @ 
PFOOOC OUC GOoOCI2OeG oC 
(a f | | tt 


SUOCUOOCCOO°OUCUOCe. 2c 


3 


2HDOQOOODOCGOBaeaeaaeae 
SOs ONNOCOGOCCeOD2steae 
DATO+TONNODNODNOLSEHM 
ND OSD SF MOMAHOO S ety OO ® 
Om at Ost OS ANTE NSTMOODWNE 
@®eeee#ee#eee#®#esee @® @® ® @ ¢ @ 
ra 2 om ime ah 

. \ )uul | 


AMON KOPF OMA MOCO MOAN SF SF 
COVCOVOOVCOOVUVOOOOO0O000O 


GooKnodocoodS54oecasee]=—] 
QIV9OM OODIOVINAOCOCOAH 
MAMACAIM TOMNWMN DO NOOIMNON 
MOORE OAINPTOWROCINGE OM 
®ee@86mUmemCcOmUC OCUClUc OClCcoOmhUlUch OOmUCUcM OrmhUCcOOCmhUCU OrhUC HOhUCUCUcPhOmhUCUCc OmhUCUcMOmhUC(U Fh 
OOVUOIVDOOOVBO0C OOOO 0VU00 
\ i | rt ' ' 


a2 


Coo 


EGAR SEE. Tae 


C.9C 


EIA PS Lee 


NUMBER 


= C.80CCOCCCD Cl 


OF ITERACTIONS 


MENIMTZES V + SETADS U SOQUARPESERRORE 


CUTPUT NOQ.4, 
DESIRED MATRIX 


IS 


6 RCwWS, &€ COLUMNS 


17, 


MATRIX NUMBEK 


FOS A 
god ond pd gard ood 
(iqarde 
Ooceace 
LA CLA AIF © 
0 WO fF © 
OODdADOM 
Aiad $f Met 
eee#e#8e@¢@ 
OOQdOO6M 
4 


LLL SF ee OO 
tt md mt CC) om 
pid } 
CeCaoasae. 
Vol ease Okan kts 
DBMOMAC 
Cregg ® 
OD MO a DO ea 
@®eoee#? @# 
COCO0OO 
1 | 


Ctrrnrwn 
md od ed Cl) et oe 
i444 } 4 
in EN sat 
MOO OO0I 
foNES atv On G92 0,700) 
MOONIA & 
CO wd ad med ed CEO) 
e®*e?e#se@®@ 
OOoCae@oO 
i { 


LUN Pema SUNS 
ed md CL, peed ed ed 
| 4 i tt 
oOoeeaS 
EX QPOW MIN 
LN TOD OUN EO 
TMmOOOD 
Mm OLA IOS 
eeeee® 
OIOCOOWS 
i 


Wi Or SEC 
ed ed ed nd 
i ee 
Sood 
MQM Of wo 
DMDOO THN 
NO OA F © 
od md OO) nd Cd 
eeee8 @ 
i ( 


MILLAN OLA 
OC) ad St mt od ot 

meets 
oqo Se 
COc$gWIM OUT 
OCMMAI@M SO 
IMME DO 
Tam COMA 
e088 8@ @ 
VOQ0000 

wid t 


Q0.800000D 00 
CRITERICN NUMBER 
CRITERION NUMBER 


BETA 


93 


dee dette del BOD? 2 
Je a1 IDS 25 


1 EQUALS 


EQUALS 


Z 
E 


4.8378730-24 
Feu 76014022 


EQUALS 


CRITERICN NUMBER 


EQUALS 


4 


CRITERICN NUMBER 
ELAPSEC TIME 


C.99 


= iee 6 CC 


TCTAL ELAPSEC TIME 


MATRIX P FOLLOWS 


6 COLUMNS 


6 ROWS, 


ly 


MATRIX NUMBER 


19a) 2 


o) 
f) 
ry 
0 
2) 
A 


nd 


CCooGe 


DCEO |G ia) 


COC eG. 


OCOGGe 


20a 12 


Ge CAC eS 


GO Cea eG 


= 


aa 
ae 


i 
Cc 
—) 


‘ 


(BY GY Eg Giles Ga i ee 


a Gi GA GE 


OTmO 


« LSOCR EZ 


0.0 


O.0 
n 


Q2O'o 


CCC Basil 


2 
Q 
0 
0 
0 
0 


COOO0 © 


MATRIX GE ECERGRS 


6 RCWS,24 COLUMNS 


2 


NUMBER 


MATRIX 


NNN NEN NNN NN BONAIO 
COOMECOCOCCECCCOOCOO 


Ceececeec aoa oqaoeoonem 
Coceecoaceoor Ccooccc 
NPTOHOABOCTMOOTOOCOOYW 
DO SIN NI COM MOP (OM O Net st Ae 
eooetreoerevreeveer eee eee oe 
COCOOOCCeGaCOCCOOOCOO 
| ! ! l | 


OAINANAIAL NRA NAINIAIOIN AA 
COOCCOCe saeco CcoCccooocs 
ococee., Ceoecoceooalbkeg 
Oocnceve  OoocCCCCOeCc © 
OOF DOW COCTOCCOHOON 
ODOOMS ONE KMNMHOOO 


NMNARMOOCHMOMODMNOONM 


eeteteteee#e8 @€¢¢ee¢ ¢ ¢@ 80 @ @ 
AO See ee 
; | | ) 


ROAINIAQUNANANIN RIAN WCU 
CeC@eCeeGeeondc CoOegcee: 
OOOO eSSeeaeace Aaaoagce 
OOCeSeCeeeacKoo Cocco 
MOWNDOCAHOMOEAO ABDOCOCO 
TNMOMYNOCeR Cm MOUONWM™ 


Ne NSE MOOS R ON OGD AHH 
eeeee#e#eeeee?eeee¢e?ee @ 
COC Te 

| ie ( 


MAUI AIQUNIO NI QUIAICI QU NI HAO eH QUAI CY 
Oocceqg@geee Coc OCcoece@ 


COCCSeS ea aeceec ec amt Ooe 
Ono. @@Ge ee mwa CCOCcCeGecG 
Oc: oOCoCconODorroo0c°ncme 
Rat POONS TON MAOAMNPE oO 
DM ARMA DAI DAIQLAAINID FL 
@®eee?eseese#ee#e#s#eee*# e# @® @ @ @ 6 
SOO ee Ce CO ee OG 
a | 1 


NOAIMOAIQEOAINQINIAL NANNY 
QOOOCOG@2ZeCce CCOOGeCcoe 
CCOCe2aeece 2c eeaa® & 
COOCCeeee oO DQlOC Oe © 
COOCODOGWMOO9O DCOCOMMOO 
DHMMMOARMDA ACAIMMAOO 


ooeoeoeo eee eee we eo eo eo oO oO 
OoCceoeeqec oOo COcooeceCco 
| i 4 l 


CUOMIONIAIOI ALAS OU NONI AL AI OI OAS et OG 
OOOO CO UOC oOCoCcooeeG 


OoniaSse ooo ooogwock 
VO VO VOVOOOCTCOVONOCO 
OQOOOmMOOCcCCC COG OO 
MOWAT ASAE OMAN OS 
AOR ODOM a mt FO OD MALAI SO 
e@eeteteee#e#e##e#e%eee## @#e¢0¢@6hucacei;ehrehlhO 
OOCO00O O@SeeOUDUGO Ce OG 
1 | ! ! l | 


NAINANNNNN NNN GONAIOA 
YB JE ing EB: 1 SS 1 a | 5 | Sg 8 I a 


SCCD ee 


OCeaaceeae&eanes]e2aaga Go) 
OOOO O Oe 2a TODO OeO® 
Ot CORDNOOODOCOTOOVOCO 
NICO ALD O LAM DANS SHO ODRON AOD 
DEN FAIMAIM FOS UU DO MOAIC 
eeeeeeee#eeeee?e#e#ee%se#ee¢@e 

| l ‘ i \ 


NIN OReANIOIAIOIONIN AI OI A SIA OI AN 
OCOVDOOVODOVCOOWCOVO00 


Oooooeo be sooaf sae e2 es) 
OOVOOOOO VOVICOCOO0 
CSCVUSKrRrOOCOCOOCOONMOO 
MmNOODTNMATTUMOTORMD 
DODA AD FIND HH OAMOAIN PIN O SF 
@®eeeee#eeeeeeee#e~*ee#e#eee#e ® 
OCOOOOCOCOCOOCee eo 
y t ' I | ' l 


94 


MATRIX €2 SE CREENWs 


6 RCWS,24 CCLUMNS 


34 


MATRIX NUMBER 


MO MO BAIS MOON MOG OH NOAF MO OES 
OOCODOOCOCOOCANSCCOOM 


Cee eSeQacooonanaoonakee 
TSTCOMDMOCDOOMC OMOETOCW 
Aaah OPF MOONE ONOSTC OOS LO 
DOOM S MOMMMEHCORNUO A 
PFO SO Set PTHIHOrmnrOn aeastinn 
@®e%&®e#eseseeesesee ¢ @ @ @ @ 6 @ 
ee IS OS Sra 

{ ( 


KHOA CAVE ANAM MOMMO MS OS 
CCGOOC COOMCCcCocoocoe 
SEeaceoao cocecececEecece 
COorTFoo OCCoOoKvOCTUeas 
CNOCHORG COONS FESOMM 
MRAM EOR CODMCOOOOFT OO 


DO AUD EEN SG GB ett BR me et LOA 4 

eoscoevoevoereevt eevee eee & @ 

COCoccaococooocweagy] 
, ' 1 t 


CODE COS SO MOAAIM HON OAOSMOS 
CCOQOCOCCCCCO Gece 
OOoOoCOoOooOOoeCEOS CGeeadge 
MONMNOOANCOOCOC OMO0e 
CMAROMNNAONIAGA OdDNMA 
HOMO ACMOGaAHGmO NOSTOO 


NS met ed ed etd MIO ON OF OO MOD SS 
eoecuoeveereee eevee ee eo @ 
COCOOCCCOOCOCOSGG esas 

{ i , 4 | { 


CA COCA S$ KK CONG SFM OQVIAKAIS 
AQGDOOCCCCOCOCOCeGGage 


CCCOHCOoC CB OOo eese = 
COFFS ORCCCHCCOe OC ae.) 
AN AFNOCCONOCOCTOUWOM 
RM OMFCRNAOMMNODANANW 
OO mt wt OP OO wet OIA SE OO A eet ed ed 
eoeetee*eeee@eee8ee%8 @ @ @ @ ®@ 
COOCCC CC? Ce oC OOeGeeae 
{ ) | ' | { 


MMA EO RMALMOAL MOAMMAMEMS 


DOCOOMOOONOC doo Wweqoqqw. 
EeaSclcocacnl ocetkea@ea a 
ODCOUaQaoCCeO Oc OOeegGmecd 
ANIDHOSTROND NTLODDNONM 
TU DMMMAINGI WOrWUCON KO 


MSFOHAIOA AIM ONCN TNA 

ee *%e*esxe%*e%*e 08 @ © @ © © @ ® 

COTOCCOCMOOCOOUCeOOOO 
{ ( | 


CS HNCVST FH AIMA MCE MOMOYMON SL 
CO OoocvvTco cnococoecr. 


Oo AcCCcoonanaooooe@eee 
Oo sTOCoerc COCOoOcCOoOeCVecg 
a TOSBHtACOICcaernsOOMN 
Mm SOR Hl NNO CMOS AGP 
mC NO OD eed ed EL eed mt OU ed nf OD et ed 
ee @ @@ @e@ @ & © @ @ eh OemhUcOhUC OhUCUC OhUCUcHhlUhO 
COOCWOVCVOCCUMOOCOGeea] 
! ' | | { { ' 


eatvataataataal salantaalantaataaiotlaalaal Teataal sa 
COOCCOCCL VUOUGOC GIES 


EY af i Yn Ya et en | Df cf a J) SS cl | 
OODPTINNOCOCTCOOOCOSBNO 
DUN FOTOKCINOCNODNOSH MK 
NAO ACT OM BOD TT FT aS O OO 
SM TOS OMAN STN OO OWF- 
e*eoeeeeeetee#e8 @®@6Uh}?OOWmhUCcChmhUCUCOCmhUCUchOmhUCUcMHmhUCUCUcMhHMhUCUCMHhUCUCOhU! 
ith « ie oe 
{ { { i 


NAKA COAST OAM AH OAC OM AION ST Sr 
COOCOVOVCCOODVOOODO0O00O 


Oooo oansooeaoosedks. 
WOVUOFCLOOCUGDOVO0ADMM! 
MOONM TOMMNOENGCMOANON 
MOMS Say GOT OCMNNOM 
DBAS mt et D AQIS cnt tO et 0 eet et ed 
e®oeete ef @ @6mhmUhOMmUCcOOhmhUCOrhUC mhUCUchH HhUhHHhUCUCAOhUhHmhUh Hh! 
riiritee 1 ma 
{ ! | eI { 


= C .U10Ge 
oO OCG 


ERRCR IN BLANCA 
ERRCR IN @C@LANCA 


9 


wn 


MATRIX U FOLLOWS 


4y 6 RGWS,24 COLUMNS 


MATRIX NUMBER 


NANA HAR ANAR A NR AAAUR 
CCOCCOCOCCOCOCCACOC O¢ 


Cees oOocCceOoOocoocecaeo 
MOA DPE WV 0 I FOS OR et rt 
NAM AIM FINE HONDO NE NE 
OOOO Mm SHINO OS SONOS 
ILA FINE OOO OE oN SON 
®e6e¢:«e8#88 @ ¢ @ ¢ @ @ @ @ « 8@ @ @ 
fae oe OO OSC OCCS 
{ { ' 


ONNANA ON BIAIOAIOI AIO AIQI OU AG 
Se] OOOCOCOCOCOOOOCSo 


Saoaeeaoeaceeaoce Cc G 
SR-OMHOTODMDEYSMODNOHO 
HCO OFTDMODOOCH ANE Hc 
OREAIOMNCMHA ARE Oana m aoe 
M™ AICO AI AO Net SOM eA C ON 
ee «© #@e#8eees*eeseseeeseseeeee ee ® 
mae Pee ee SO CO 
{ j i! 


NOUNNIALAIOQI A et QUAI OI COON ot QUAL 
Geeceocc OC CcCoeocooecCc 


OoOeee:Cc$ockooo@Geqgeacea 
aot COME OO OE CMOORAMNN 
CO mt =P FR LINAC OOF COOP 
Det SOHO AC ANAIS OO OMMMINA© 
SaaS ™ OO SIN OS ep OO Meaty 
®e#etee#e#ee6¢#208 84 8080800800 «@ 
Seo CCOOCCOCOOCOeCcC eG 

I { {tet { 


NANA ONAINAIN AeA EAA AN 
COococecocoocooeocCcoce 


SSeS jOOoceaeageeoga] 
Forts Morn tem ON—aeoO OP 
at OS USO OSINSHIOSE COO 
CA COM LAM DA SOU OD HO & ot 
DOHC M$FOCST AO AI SAARI LA SO 
e® e0@«©@ © @6mhUhUhOmhC<C MHOC TC OlhC<C HTC K HOChUCUhOOmhUCU OrhC OrmhUCHDhC—<C NM HhC<“C MCSFHhClUO 
ip) Vises 
| { { 


NNN AIOININ AIRS AIAINN OA a 
COCCCOOOOCOCONAICOCUCS 


{ 
Seen CeeolCOCaCeCaoCac 
O NOM OFTGO OF MOMANN OW 
ORADONADOMOKOMANHt OS 
POD SIN A LD OO OA OOM 
eet @ G8 8e@8mUMOMmUMOMWUCUOHMUCUCOHMUCUCOUCOmUCMHMhUhHhUhUh Oh 
(ey Fd J) ot RT <a ce I ee GD I eae Fay I ib] ean I 2] a Yi ato 7 aa. 0 
{ { | { 


NO NINOAIN BAI AIOAQINI QUIN ANIA AN 
COODOVOOCOONOCOCCVOCOC 


SSS OoaaQahoooncoeceaG 
DORAIVA WAT AN CR OD aA M 
WO MP AININ OA ©O O TLE DINO AD 
OCOMaABOOSAMOAM-PES MOOD 
NOUN HO COND OH) HONG OMY AI SO 
® e0e@B ® @@8eeeee¢6.e8thmhUhthmhC OhUC MUCrh!hUrhhUhFO hO 
OGeoCOOCCOC OC eoOoOeoCce eG 
{ { | { : { | 


NAINA AIAQINIOQUN NA eSOAIN ON ANA 
COCODSGCCCOSDOVOOCCC 


SOoSQaogaqggQgaQqgtOooagageagm 
DMrATOODVOSNOM That or 
INO DO BOND ALO OAD QD Ht DUN FP O 
0 OHO DOMINAN SHO SAM et OD 
DRA S NOAA SF OS NAM OMAIN 
esooevreeoeene ee eee © & 8 
ae te ae Oe OOO 

| { { pet 4 { 


MAIOMRAIAIAIONAIN AIAN Ra ANAIAIO 
OOCOOVOOCOCOOCOO9O00 


A sf if Yes fy I ea Tn Fa ef a Yate Fc 

NTO TONRAeHNOOMMOPROTN 

US CAIM SRANIM AIM MMINA SAT in eS 

OP HARIBO OOOMOSOCNDAMOW 

MAINE OINMING HORA SIN OS 
®eeeeee+eeseeees#eseee#?#s 8@® @® ® 

OMOCOOCCOCVVOCOOOO00O 
| { { '4 \ 1 


MATRIX V FOLLChS 


on é RCwWS,24 CCLUMNS 


MATRIX NUMBER 


MAAC A TMOOAMNMAM MO CAIM SOO 
OOOSCOOOCONGCCOCVOOCOO 


Seaocococooocoooet ees. es] 
BOOKMAILNAMNOOMHOMOO 
MNOSTASNOMNAHOOAMNM ONAN 
SHON CM HMALPNDOAKNN aS 
FINOD Smt Feat OP MOON ett OS 
eeoerevreeweeeeeeve eee @ @ 

' | 


FAAINMNOMONNMAMOAMA MOOS OM 
OTC OTCOOVOODVDVOOOGO0COO 


Saccooocnaoacneeeseeae]} 
BASHAM TINGOLTPROHMNAMNS SM 
MADONNA OHO GS OW 
HINA COMmORrOHOaAOOONS™ 
DO AIC COD LA OF OO wt oO ed ett OU OF 
@e*eteee@eee%ke#eee?@¢¢e0e80 @ @ 
OSL Ole! e111 Ole] Ole ®lelelelelelele) 
| | 1 1] 


AO MOOS FOAMOAANMOAMHRMOMAIST MH 
ODOSCSCMOCOOOCCOOCOCOoOe 


OCcCeSconaoeeeeaceewewe Ss 
ar NSCHAODSTOMOCOMNG!> 
CAF OSHA OM ADAMO Pr ONO 
MOOP SCOP MH AW OW Nf OF 4 
ONY et ON eed red ad md CI LD SO ed OD ee OS 
eetetee¢e¢@ee8ee¢eeee¢ @¢¢@ @ @ 
SO OS OO a aaa ie 

f | | 


AMA SF FAIOMOAMN SF OMAINA St 
DO COCMVWDWDOONOVO00°O 


Ce SCOoonocOoecaoeeaea 
OM NO-EMANNCSOMAMDOSEH AN 
CHAM ONAN SONNAMMODMONST 
NOAA AOMMAAININO GMOS Ct~A4A0 
DS aA et OM AS HAIN OS RH 
e@eeese#eee-eseee#e@eee0ee0e0@8086h%3 @ @ @ @ 
CSOCOCO OOO CCOCeeeoOe 
l 1 | 1 | | 


PAO NP OFANANOMNAM OOS MO 
OOTOOOVGOCOOCVMVOCOOOO 


CCococncooneeacke aka se] 
PANTOAMN CHAIN OCODAN Om 
TMH AMOOMNDONDADORHOMC DSO 
COON OAMOMTOONALIEMNOC 
OSTA OAM ONE UM AIAN = 
eoeteee¢e¢eteeeeeseeee¢e @ @ 
Coco CoOOo GoOocooooeGd. 
I it ‘ 


CVD OVCE SAP COON CRMC On eee Gain) aa 
OUDVVORBVOOCCOVUGOODVLO00 


TNO OMT ADO ANIAN DO OANA TC 

BOHM POOODOMT DOCAYQO 

WMI AIA C4 OO GoM NO Ooo 

4) DO Stet PTE PANS 
eeeeereee*e@eeeeeeeee @ e 

COC OCCOOOCO BOC eoeeeo”m 
| i | | ! 1 , | 


COI OT ONLY CES. COON OES OV Cree ta) oe 
OVUCOSC OOVOODCOOVOOOUO 


OHOOodoOo0o0OOBoCeaSeeaSss 
NOD PLN I POA PTR MAT HAHOO 
SAMY OOWUOMN OSIAMOMONID 
OMIA NOAA AIA HANLIS OM Eh 
SO PM TO Sete SINAN OOM LAD 
*eee#eeeeeeee20e¢0¢06h6Uc3ca @80 0 060 086 «= 
Oe Saale | a a 
| | ‘ 1 


OAM AOS TOMA OHOACOCYIAIIO SF SS 


OoOoCooO CoG oaeaOooegeae 
QOT AO “AUVUID OM WU NUN STP 
ODM OO TINA GAD AOM MRA LAN 
ARONA DOOMOHNT Deh NNOAI 
AALS MH aA NQ\N OAH AN AMO Ga 
@eeeeeseeeeet¢?#ee% 8 © @e @ 
YQOOVOCODOVOVIOOVDO9D00O 
| | | i | \ | 


96 


FLAFSED YRME 


1.66 


ELAPSED TIME 
NUMBER OF ITTERACTIOQNS 


Gel7CCOCCcOD C2 


MINIMEBZES V @ BET& * U SQUARESERRDRS 


CUTPLT NG.4, 
DESERED MATRIX 


OS 


6 ROWS, 6 COLUMNS 


l7, 


MATRIX NUMBER 


COCOAN 
COOOCOO 


{ 
COQAQ0GQAG 
TOAD eth a 
onto ok loa r of me 
MAM MOO 
aaOANIOAl 
e088 @ @ @ 
ar ea 

| { 


COOAHOxt= 
ONdO00OO 


! ! 
CoeeacGce 
INNO CUS 
LAS NOM 
LAD FOWM 
aA OS 


OOONoOo 
OOCONOO 


CGoeaasoO 
ONACOs 
UVa R= Oey 
SNOMHOr 
MMO ANA 
eo ¢?0e080 @ ® 
OdQTo00OO 

\ 1s 


aAOQOHOHO 
OQOUDDO 
| | 

CaO Grec 
ONADOW 
“= LAO AIO 
m™oOooeoysre 
tot TOUR 
e*¢37oe@¢ @ 
COOVVMOY 
| ( 


ONOOMO 
OOO OOO 


Ores ee) 
DHraniayanyad 
TO TY OV LN SISO 
MP OO ON 
md ed md CO) end med 
eoe3ee8e?8 @ 
OVOOOoO 
| | | 


aOAOOO 
OVOOVO 


' 

SO ar) Oe 
TAoownr 
OW COUN 
NOM AQUA 
OD afd CO) nd md 
oe 80h68hUmhlUO 
Oaceoo 

4 | 


O.8900C0D O00 


BETA 
—_CRITERICN NUMBER 


4.S98249D 05 


1 EQUALS 
2 EQUALS 
3 EQUALS 
4 EQUALS 


97 


2-57C142D 03 


“a 


CRITERION NUMBER 
CRITERICN NUMBER 
CRITERION NUMBER 


139593D 94 
22188228) O05 


a & 


C0.9S 


— 
— 


ELAPSED TIME 


ieee ts, 


TOTAL ELAPSEC TIME 


APPENDIX > 
BRIEF TREATMENT OF RELATED PROBLEMS 
D-l1 Case where there are no Restrictions on C. 
In the development of this thesis, starting with the matrix equa- 
tion CU = V, matrix C was restricted to be symmetric and positive definite. 
Suppose now that these restrictions should be removed, and, indeed, let 
it only be required that C have the correct dimensions to relate matrices 
U and V which may have a different number of rows. 


Let us write equation (4) in a modified form 
T 
A=t[(cu-v) Z(cu-v) | (D-1) 


where Z is arbitrary and might be thought of as representing a set of 


weighting factors. If Z is the unit matrix, equation (D-1) reduces to 
equation (4). 


From equation (D-1) 
K=h [UC ZCU-U'C'2V-V'ZCU +V ZV | 


Note that, in general, cha C since the requirement of symmetry 


is not made here. 


Using * to denote partial differentiation with respect to any element 


x, as before, one gets 
X*2 te [UTC *ZCU UC 2ZO°U -U'C *ZV-VZC*U | 


=th [CUE = Ve VR eaieg hea eun | (D-2) 
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Thus CA is made stationary with respect to cach element of C if C 


satisfies the equation 


If U is of rank m, where m is the number of rows in U, then 
Coe vu (uuty 7} 
It is interesting to note, that the result is independent of the 


arbitrary weighting matrix Z. 


D-2 Case where C is Orthogonal. 
This problem arises in the treatment of the orientation of a rigid 
body. 


If C is orthogonal, then 


~1 « e e Ld e 
where the notation Y = C is used for convenience in manipulation. 





As before let x be any element of C, and let * denote 2 
Then the relation CY = I, leads to C* = -CY*C, and equation (4) can be 
written as 
A= tr (UU-U'YV-V'CU +V'V) 
so that 


Ae — lial OY Y + Vey =aU } 


=th[(Y")(CUv"c -vu")] (D-3) 


Thus, & is stationary with respect to variation in any of the ele- 
ments of C, if C satisfies the equation 
cuv'c = ww! 


which will be so if C is any of the several values 


c = (vutuv’)*(uvt) 7} 


oo 


vutuv? is symmetric matrix, and its square root can be obtained as 
described in Section 4. 
These results were obtained by Professor Brock and are given here 


so as to collect in one place information on this general matrix problem. 
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